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ABSTRACT 


A digital computer program is presented capable 
of determining the imand/or out-of-plane vibration 
frequencies of a single plane piping system. The 
analysis is made utilizing the method of transfer matrices, 
a discussion of which is included. The program accepts 
branched systems with or without intermediate supports. 
The effects of both shear deflection and rotational 
inertia are considered although either or both may be 
neglected at the user's discretion. Distributed mass is 
generally treated as such, although provision is made 
for alternate treatment by lumping mass, 

Several typical viping systems are analyzed using 
these various program options, and the results obtained 


are compared. 
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CHAPE A Rel 


INTRODUCTION 

1.1 General Remarks 

The dynamic behavior of piping systems, particularly 
their response to periodic and/or transient excitation, 
ils a subject of considerable interest and growing 
pertinence in the design of such systems. The reasons 
for this are readily discernable when the unprecedented 
reliability demanded of nuclear power plant piping, or 
the consequence of a piping system failure in an orbiting 
space vehicle, are considered. An increasing number of 
piping systems are being designed, of which, like those 
mentioned above, a long service life is demanded, and 
to which access is limited or completely denied once the 
system is placed in service. Clearly, no factor which 
might contribute to the failure of such a system can be 
ignored. Fatigue brought about by periodic excitation 
at or near a system resonance is one such factor. Cata-= 
Strophic failure caused by a periodic disturbance of 
Short duration and high amplitude, such as an earthquake, 
is another. 

With regard to the latter case, it has come to the 
writer's attention since beginning this work, that 


Japanese engineers are giving considerable and close 








attention to this particular problem. Their first step 
in approaching it is to determine the natural frequencies 
of the system of interest. 

Fortunately, or perhaps fostered by necessity, 
development of the “tools" needed to carry out accurate, 
reliable, and practicable vibration analyses of complex 
systems, namely analysis by matrix methods with the aid 
of high speed digital computers, has kept pace with the 


increasing demand for their employment. 


1.2 Scope of Work Presented 

A method for determining, analytically, the in-and 
out-of-plane undamped natural frequencies of planar 
piping systems has been developed and is presented 
herein. It employs the method of "Transfer Matrices"! 
to construct an accurate mathematical system model and 
determine system natural frequencies, and a digital 
computer to carry out the inherent voluminous numerical 
computations. VIPIPE, the digital computer program pre- 
sented in appendices A through C constitutes the embodi- 
ment of the method. 





ithe theory and method of Transfer Matrices is dis- 
cussed in sufficient depth for the purposes of this 
report in chapter 2. A more comprehensive and detailed 
treatment of the method may be found in ref. i. 








This method is a first but important step toward 
the ultimate objective of analysis of three dimensional 
systems of unrestricted topology. 

Attention has been confined to the two dimensional 
or planar case because the theory underlying the treat- 
ment of three dimensional systems using Transfer Matrices 
1s not presently available in the literature, and because 
it is felt that the need for, and utility of, a treat- 
ment of the two dimensional case fully warrants the 
development. 

Although in its present state of development, the 
theory of Transfer Matrices permits the consideration 
of any planar piping configuration, certain restrictions 
are imposed on the configurations which may be treated 
by the method herein presented, by the limitations of the 
program VIPIPE. These limitations are discussed in detail 
in section 1.3. 

The philosophy of VIPIPE has sufficient generality 
to make the program readily adaptable to consideration 
of the three dimensional case as soon as the underlying 


theory is developed. 


1.3 Assumptions and Limitations 


The mathematical models used to describe the 
elastic properties of a system and its components are 
faithful within the framework of the following con- 


Siderations : 








a. The material of which each conponent is made 

is considered to be homogeneous, isotropic, and 
linearly elastic. 

b. The system is conservative (no damping). 

ec. Distributed mass is treated as such except in 
curved sections for which a lumped mass model is 
used. For curved sections, sufficient sub- 
sectioning is carried out so as to reduce error 

to an economicallytacceptable level. 

dad. Spring hangers are considered to consist of 

from one to three mutually perpendicular linear 
springs, and from one to three torsional springs 
which act about mutually perpendicular axes. When 

a@ spring hanger consists of both linear and torsional 
springs, the linear springs are normally considered 
to act along the same mutually perpendicular axes 
that the torsional spring act about. <A more gen- 
eral case may be treated with relative ease and 

is a subject of discussion in appendix B, section B.960 


Although the theory would permit determining any 


number of natural frequencies for any planar system 


having any number of components, branches, and hangers, 


program VIPIPE, described herein, is limited to the 


treatment of systems having a maximum of fifty components, 


lappendix B, section B.93. 








twenty branches, and twenty hangers. Further, only 


Single filament branches from a main member may be 


treated, i.e., all branches must emanate from the same 


continuous member, and a branch may not emanate from 


a branch. 


1.4 Netation 


Lal 
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u,V,W 


Ry fo 


op) 


matrix 


‘cartesian right-handed coordinate system 


displacements in the x,y,z directions 
respectively 


rotations about the x,y,Z axes respectively 
displacement vector 

Young’s modulus 

pending stiffness 

Shear modulus 

torsional stiffness 

coordinate transformation matrix 
unit matrix 

polar second moment of area 
length 

bending moment 

normal force 

null matrix 

internal force vector 

square submatrices 


spring matrix 








torque 

transfer matrix 
transfer matrix element 
shear sionce 

unknown vector 

state vector 


subscripts, left and right side of com- 
ponents respectively 








CHAPTER If 


TRANSFER MATRICES 


2.1 The State Vector and Transfer Matrix 
a. State Vector The state vector at a point i of 
an elastic system is a column matrix having sufficient 
elements to completely describe the instantaneous dis- 
placement and/or change in orientation of that point 
from its quiescent position, as well as the force and 
moment system existing in the member at the same point 
at the same instant, i.e., those forces and moments 
which, if a cutting plane were passed through the point 
at the instant of interest, would have to be applied to 
each cut face to prevent relative motion between then. 
For perfectly general structural systems the 
notion of a state vector is a rather sophisticated 
concept. In any properly constituted state vector 
however the number of elements pertaining to forces and 
moments is equal to the number pertaining to displacement. 
In the case of a planar piping system, the state 
vector has six elements, three pertaining to displace- 
ment, and three pertaining to force. The state vectors 


i 


for the im and out-of-plane cases“respectively, as used 


luenceforth the in-plane case will be referred to 
as Case 1, the out-of-plane case as Case 2. 








in VIPIPEK, were chosen to be 


aa 67 
Vv < 
Ty -W 
Zip= My Zoop= | (2.1-1) 
Vy My 
N V 


The elements of these state vectors are the magnitudes 
of vector quantities of the same name which derive their 
orientations from those of tae X,y,2 axes of aliiocall 
right-handed cartesian coordinate system (fig. 2.1-1), 


associated with the particular point of interest. 





Fig. 2.1-1 Coordinate system and vectorial representa- 
tion of state vector elements. (Interpret double shaft 


arrows by the right-hand rule.) 


These local coordinate systems are constructed 
such that the plane of the quiescent system and the 
xy plane of the coordinate system coincide. The local 
xX axis is coincident with the local tangent to the 
quiescent system centroidal axis. No reference is made 
toa common fixed coordinate system. In establishing the 
coordinate system for the purpose of employing VIPIPL, 
the longest continuous run of pipe, or the one from 
which all branches emanate if they are not one and the 
same, is designated as the main member, one end of which 
is chosen as the starting point. This end is called 
the left end and its state vector is symbolized Z, . 
The positive x direction is taken as away from the "left" 
end and toward the "right" end along the centroidal 
axis of the main member. The positive z direction is 
perpendicular to the quiescent plane of the system and 
invariant, and this direction is called downward. The 
y axis is chosen so as to complete a right-handed car- 
tesian coordinate system. 
b. Transfer Matrix Let us consider the relationship 
between the state vectors at the left and right ends of 
& massless straight section of pipe. We begin by con- 
sidering the Situation occurring in Case 2 in which 
both bending and twisting (about the centroidal axis) 


occur. The relationships between the left and right 








end displacements are given by; 


= lgle e 
&. = S sacle (2.1=2) 
- = + WL +M ia + V pica (2.1=3) 
We = Wp 4 ET ere 
_ eel (2.1-4) 
URE Mie LES OL ¢E 
The force relationships are given by: 
TR = Ty < (2.1=5) 
Mp = iy = V,L (2.1=6) 
VR = Vi. (2.1=-7) 


or, in matrix notation 


fe eK 0 ; 
GJt 

rc Od OO CnO 7 
-w o o 1 L BS LF lew 

= 2EJ 65 
alg a es ey 

“he thd Y 
My O- 0 “0 “*OeeieaaE My 
V, ad) 0 GOnom dt Vv 
Z R - Z L 


The six-by-six matrix relating the state vectors of 

the left and right end of the pipe is termed a transfer 
matrix. It is the transfer matrix of the elastomechical 
system constituted by the pipe and its boundary condi- 
tions. Similar sixeby-six transfer matrices may be 


developed for point masses, curved sections of pipe, 


10 








and elbows, etc. Many of these are contained in the 
transfer matrix catalogemicluded in ref. ii and are 
used in VIPIPE. 

Since the application of interest is to the deter- 
mination of the natural frequencies of piping systems, 
it would be instructive at this voint to illustrate the 
origin and nature of frequency dependence in transfer 
matrices. To that end, let us construct the transfer 
gatrix of a point mass, again considering Case 2. To 
begin, note that throughout the application considered 
herein we are concerned with the isochronous vibration 
of a system at a single frequencya. The tacit assump- 
tion has been made that all quantities varying with 
time, i.e., all of the elements of the state vectors, 
do so with the form; constant « coswt, and have exhibited 
only the constant or amplitude, while suppressing the 
variable factor. 

The transfer matrix of a point mass may be constructed 
by noting that deflection, twist, slope, and moment are 
constant and continuous across the zero length section 
which a point mass constitutes. The vibrating mass does 
introduce a discontinuity in shear however, which may be 


determined by considering the free body diagram of 


1See App. F. 
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Fig. 2.1-2, and the following relationships; 


w(t) = wecoswt 


2 Bs 
a ae 
we) WW COSmt (29) 
at 
From F=ma; 

c 

V, Roost -V,, cosat=and a (2 a=d\0) 
dt 


and substituting (2.1-9) into (2.1=10) we have; 
Veg = Vez + mo*(-w) (2.1-11) 


and are ready to construct the transfer matrix UL of 


a point mass. Hence 


1 0 0 0 0 0 
0 1 0 0 0 0 
Up = 0 0 i 0 0 0 
(2.1-12) 
0 0 0 ut 0 0 
0 0 0 0 1 0 
0 0 eS 0 0 1 
im 
sale View 
mos 


Fig. 2.1-2 Free body diagram of a point mass 
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Consider now an elastomechanical system, such as 
that shown in Fig. 2.1-3, made up of many components 
connected end-to-end, each of which, as should now be 
apparent, may be characterized by a unique transfer 


matrix. Subscripting the 


Fig. 2.1=3 An elastomechanical system 


transfer matrices to reflect the subscripts assigned 
to the state vectors which describe the force systems 
and displacements at the ends of each component, we may 


write the following matrix equations; 
42-04 2°41; 43=U2 3°22 50005 4n=Unn-1 n° 4n-1 (2.1-13) 


Note that successive state vectors may be eliminated, 
by substituting for each its matrix product equivalent, . 


until finally the matrix relationship 
4, =Un-1,n°°°U2,3°U1 2°24 remains. The matrix product, 


Un-1,n°°°U3 2°02 4 , of the component transfer matrices 
is itself an m-by-m matrix which relates the final to 


the initial state vector of the system. It is, in fact, 


the transfer matrix of the system. 


NS) 








2.2 Hlastic Intermediate Conditions and Branched Systems 


Piping systems frequently include branches, and 
are supported intermittently along their length by hangers. 
These branches and hangers usually have a very signifi- 
cant effect on the dynamic behavior of a system. Their 
effects may be considered within the framework of the 
transfer matrix method. 

a. Hangers It is clear that the effect of a hanger on 
a system is to introduce discontinuities in the internal 
forces in the system at the point of attachment. 

Consider a spring hanger having the torsional spring 
constants CTX about the x axis, and CTY about the y axis, 
and the linear spring constant CLX in the x direction. 
For a Case 2 situation, the equations relating the state 


vector elements on either side of the point of attachment 


ares 
feiss 
TR = Ty + CTX*°s 
-Wp ==Wy) 
We = Wy, 


<j 
wy 
M 
ue 
+ 
a 
cS 
o 
= 
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or in matrix notation; 


is 1 0 0 0 0 1s | 
T cTx 1 0 0 0 T 
=-W 9) U 1 0 0 =W 
= |0 0 ob at 0 y 
i 0. 80> WG Remy oer ol | M 
V. 0 0 <=CLX OoOMee 1 y 


It should be apparent at this point that a rigid 
intermediate support, i.e., one having one or more 
infinitely stiff springs, may not be considered using 
this method. Rigid supports may be approximated however 
by uSing very large values for appropriate spring con- 
stants. A value of 100,000 lbs./in. or in.=-1b./ rad. 
as applicable would be an adequate approximation in most 
instances, and it does not introduce computational diffi- 
culties. 
be. Branches The joining of a branch to a system 
has an effect Shumieee be that of a spring hanger on the 
force system in the main member at the point of joining. 
In fact, it “looks" exactly like a spring hanger to the 
main member except, as will shortly be seen, its parameters 
depend on frequency. ODiscontinuities are introduced in 
the forces whereas the main member deflections on each 
Side of the joint are the same, and identical to those 


of the branch at the same point. 


15 








Consider a Case 1 situation and systen shown in 


Bie. Areca lee 


Lh 7 


li, ek 
ane system 
B 
[a 






Fig. 2.2-1 Main system with a branch 
In the free body diagram, Fig. 2.2=2a, of the joint of 
the system of Fig. 2.2-1, the force systems acting on 
the joint are expressed in terms of their own coordinate 


: By a suitable coordinate transformation it 


systems. 
is possible to express the force system of the end of 


the branch in terms of the main 







M a Mp Me Me 
(| mn, oy tN 
&/ Vp a; ‘R 

Nate 


Fig. 2o2-2a Fig. 202-2bd 
member coordinate system. Such a situation is shown in 
Pig. 2.2-20. 

lfhe symbol ~ refers to the fact that the vector 


@lements are expressed in terms of the coordinates of 
the branch. 
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Equilibrium considerations and the orthogonality of 
the force and shear vectors permit us to write the 
following equations: 


Y= Mee 


R L B 
V=V +V 
R L B 
N =N +N 
R L B 
or in matrix form 
M 
+ = (2.2-0) 
h N N 
L B R 


Bearing in mind the continuity of displacements 
at a joint, it is clear that the main member state 
vector Zn can be obtained by simply adding the force 
components of the branch state vector 23; expressed 
in the coordinate system of the main member, to the 
force components of the main member state vector 4° 
What follows is the derivation of a suitable expression 
for the force components of the branch state vector. 
The main member state vector a is obtained easily 
from the state vector aN and the transfer matrices of 
the intervening main member components. 

At point C (Fig. 2.2-l)there are three zero and 
ennee non zero elements in the state vector Zo» 

Isee Appendix B, section B.92 for treatment of end 


conditions which do not meet the requirement for three 
zero and three non zero elements. 
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The column matrix of the non zero components is 
symbolized ¥. after relating the branch state vector 
at B to that at C by the transfer matrix of the branch 
(obtained in the usual manner from the transfer 
matrices of its components), and eliminating those 
columns of the branch transfer matrix which would 
otherwise be multiplied by the zeros of Zos we obtain 


the matrix equation 


aa v (222=40) 


wherein the state vector Zp has been reordered and par-= 
titioned such that d, denotes @ sub vector containing 

the elements of Zp pertaining to deflection, and Pa 

Sub vector containing those pertaining to force. Sxoan- 
Sion of equation 2.2-1 gives 

do =e ea (2.2-2) 


and 


Pa = BR. ¥ (2.2-3) 


By substituting ae : dp = y in equation 2.2=3 we 


obtain 


Py = Reg er) al (2.2-4) 


4 B 
relating the displacements at B to the forces in the 


branch at B. The matrix product aig is called a 


spring matrix because of the function that it serves. 
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Although any consistant system .ill do, it is con- 
venient here to establish a particular nethod for measuring 
the angle at which a branch joins the main member. The 
one defined is consistant with that used in VIPIP£. 
Accordingly, the angle at which a branch joins the main 
member is measured counter-clockwise from the x axis of 
the main member as it is oriented immediately before 
the intersection, to the x axis of the branch as it is 
oriented at the point of intersection. The angle @ in 
Fig. 2.2-1 is drawn in conformance with this definition. 


From the matrix relationships 


dp = Gy? d, 
and 
Pig = Gees 


where for a planar piping system 


cosg -=sing 0 thea () 0 
G, = sing cosg Q and Gy = |0 cosg -sing 
0 0 1 . O sing cos¢g 


the spring matrix (symbolized S) ‘ 


1, G 


S = G.° Ro° Ry 1 
is obtained which relates the force and displacement 
elements of the branch transfer matrix 25 (expressed in 


terms of the coordinates of the main member). Thus 


Pp = S° dp 


deg 





It is now a simple matter to construct the transfer 
matrix relating the main member state vectors on either 
side of the joint. Recall from equation 2.2-0 that 

PR = Pie 2be= ieee 
since 


dp = dy = dp) 


' 0 ™ rol or Za = Upe ay, 
45 = pote as 
t 


t/a 7 


where we use the symbol Up to denote the point matrix 
which has just been derived. 
Thus a branch is treated as a point transfer matrix 


in constructing the transfer matrix of the system, 


2.3 Mode Frequencies and tie Frequency Determinant ~ 

A principal application of the method of transfer 
matrices, and the one with which we are concerned, is 
to the vibration analysis of elastomechanical systems 
for the purpose of obtaining their natural frequencies. 

When transfer matrices are constructed for systems 
having Mass , certain of their elements are frequency 
Heston: (eeeel1 the transfer matrix of a point mass 
derived earlier). Proper anplication of abe systen 
boundary conditions to these transfer matrices yields 


a matrix of order n/2 (from a transfer matrix of order n) 


which can be used to obtain the natural frequencies of a 


20 





systez. Its determinant is termed the frequency determi- 
nant of the system. The conceot is best described by 
an example. 

Consider a case 2 situation wherein a system has 
a system transfer matrix Us and boundary conditions 
embodied in its left and right end state vectors 4, and 
Zao Assume that the left end is fixed and that the right 
end is free. The matrix equation relating the left and 
right end state vectors is 


§ ili ° ° ; Be | |e 
° ° ° o 9 Qo Oe 
-W @ ° C.) ° © C) QO 
= 8 
° ° ° 8 ° c+] | 0 
u | 
0, ° ° ° ° °o i<] M 
Olp U6,1 ° ° ° N66 hig 


from which may be constructed the matrix relationship 


e 2,2 Yes Ye ir 
0 = Ws 2 Us i5 Us 6le M (2o9-1) 
0 R U6 ,2 46,5 46 6 V., 


Clearly, a non-trivial solution requires that the third 
order submatrix of equation 2.3=-1 must vanish. The 
determinant of the submatrix is the system frequency 
determinant mentioned above. The frequencies for which 
its value is zero are the natural frequencies of the 
system described by the transfer matrix Ug and the 
boundary condition state vectors, 4, and Zp, from which 


itewes derived. 21 





2.4% The Digital Conputer and the Method of Transfer Matrices 


Even for a relatively simple system, the elements 
of the frequency determinant become unwieldly polynomials 
in radian frequency. Determination of the roots of the 
equation resulting from an expansion of the frequency 
determinant is usually quite impractical and recourse 
must be made to a numerical solution. Ina numerical 
solution the frequency determinant is constructed and 
evaluated for selected frequencies over a range of 
interest. Its value is plotted as a function of fre- 
quency and the natural frequencies of the system which 


it represents are the zeros of the plot (Fig. 2.31). 
+ 






+09 
Lis the velue & 
the Frequency detey. 


minanls 


frequency (vad. /se 


Fig .203-1 Frequency determinant as a function of frequency 
Use of a digital computer permits rapid and accu- 
rate analysis of complex systems. Program VIPIPE 
presented in Appendices A through C, for instance, will 
accept for analysis planar piping systems having up to 
fifty components including as many as twenty branches 
and spring hangers. It treats distributed mass as such 
and considers the effects of shear deflection and rota- 


tional inertia. Ina typical analysis, the first five 








natural frequencies were found for a system having 
thirteen components, including two spring hangers and 
two branches, at the cost of five minutes of machine 
time on a CDC 1604 computer. 

It is throuch the use of high speed digital 
computers, that the full potential of the method of 
transfer matrices in this application as well as in 


otherslmay be realized. 


lsee ref. 1 fora comprehensive treatment of 
applications of the method of transfer matrices to 
other problems in elastomechanics. 
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CHAP EARS Li 
DISCUSSION 


3.1 Accuracy of Method 


Confidence in the accuracy of the method was 
established by using it to determine the natural 
frequencies of several systems for which these fre- 
quencies could be obtained by other means, 1.¢€., by 
means of either a closed form analytical solution or 
by recourse to the literature (ay. The average dif- 
ference between system natural frequencies obtained 
uSing VIPIPE and those obtained by analytical means 
was found to be less than .0064%, while the average 
difference between system natural frequencies obtained 
using VIPIPE and those obtained from the literature 
was found to be less than .5%. 

Details of the systems analyzed and the results 
obtained may be found in Appendix E, section E.21. The 
analyzed systems consisted of single spans of straight 
pipe having various combinations of boundary conditions, 

The published papers of E. Voltera [4] eve 
methods for determining the natural frequencies of 
curved filaments, but the methods are sufficiently 
involved as to make their application in obtaining 


Some comparison frequencies a major undertaking. 


1yumbers in brackets pertain to bibliography. 
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Time precluded assumption of the task. The literature 
[259 does provide means for computing the fundanental 
frequency of partial rings with fixed ends. The dif- 
ference between the fundamental frequencies determined 
uSing VIPIPE and those computed using the method of 
the literature, was found to be less than.3% for in- 
plane, and 2.3% for out-of-plane vibrations.! Details 
of the system analyzed and the results obtained may be 
found in Appendix E, section E.24. 
3.2 Integrity of Solution 

Confidence in the integrity of the method was 
established by conducting"mirror image" analyses of 
almost every configuration considered. Construction 
of the transfer matrix of a system is directional in 
nature, in that an end of the system is chosen as a 
starting point, and the transfer matrix of the system 
is developed by premultiplying the incomplete system 
transfer matrix by the transfer matrix of each system 
component in the order that each component is encoune= 
tered in tracing the system from that starting point. 
It follows therefore that an error in the structure of 
any of the component transfer matrices. or in the method 
or order of their combination in constructing the 


system transfer matrix, would reveal itself in the form 


tee footnote to section E.24 of Appendix E. 


nd 
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of inconsistent results if the system were analyzed 
twice, beginning the first time from one end, the second 
time from the other end. Such a dual analysis is 
referred to herein as a "mirror image" check. 

In general, the disagreement between natural 
frequencies obtained for a system and for its mirror 
image was found to be less than .0003/. 

An added check on solution integrity was obtained 
by analyzing a single branch system treating first 
one member then another as the branch. The average 
a@ifference here was on the order of .027% for in-plane 
vibration, and .007%4 for out-of=plane vibration. 

Details of the systems analyzed and of the results 
obtained may be found in Appendix E, sections E.31 and E.32. 
3-3 Computer Imposed Limitations 

The number of natural frequencies that may be found 
for a given system using VIPIPE, is a function of the 
floating point significant figure capacity of the com- 
puter used. The significant figure capacity required 
to solve for a given number of natural frequencies of 
any system is a function of the characteristics the 
system; and for all systems, the significant figure 
capacity required increases with the number of the mode 
sought. Thus, for example, for a given system, the 
number of significant figures required in order to 


determine the third mode frequency is greater than 


26 





= - 
— 7 — 
—_ 
— eS 
a 
eee Smee — 





that required to determine the first node frequency, 
and while the first ten natural frequencies of one 
system may be obtained without difficulty, the signifi- 
cant figure capacity of the computer may be inadequate 
to determine the fourth mode frequency of another. 
Program VIPIPE is designed to terminate a solution 
when the number of significant figures required for its 
continuance exceeds that available in the machine. 

The use of multiple precision arithmetic would 


extend VIPIPE'S utility considerably. 
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APPENDIX A 
DESCRIPTION OF PROGRAM 


Aoi General Remarks 

Program VIPIPE is a CDC FORTRAN 60 [5] language 
digital computer program designed for use in the vibra- 
tion analysis of planar piping systems. In and out-of- 
plane undamoed natural frequencies may be found through 
its use. The analysis is performed by an application 
of the method of transfer matrices wherein each component 
of a system is characterized by a matrix called its 
transfer matrix. The entire system is characterized 
by the system transfer matrix obtained through an 
ordered multiplication of the component transfer 
matrices. Considerable flexibility with regard to 
the mathematical model employed to describe the systen, 
and the nature of simplifying assumptions accepted in 
its vibration analysis, is provided in the program. 
For instance, although the method of solution yielding 
the greatest accuracy coupled with least machine time 
combines a distributed mass treatment of straight 
sections with a lumped mess treatment of curved sections, 
and considers the effects of both shear deflection and 
rotational inertia, the user may elect to use a lumped 
mass treatment of all sections and to consider or neglect, 
Singly or together, the effects of shear deflection and 


rotational inertia. 
A-1 








Problem solution proceeds by an iterative process 
with radian frequency the independent variable. The 
Sign of the frequency deterainantlis used to control 
the search for and convergence to the system natural 
frequencies. The frequency increment is computed 
internally, as is the solution acceptability criterion. 
Both are accessible, however, and may be controlled by 
the user.“ Subsectioning of components treated by lumped 
mass methods is carried out internally as a function of 
component L/D (length-to-diameter) ratio and the highest 
mode sought. 

The program has the capability to analyze piping 
systems comprised of a maximum of fifty components, and 
having as many as twenty branches and twenty hangers. 

A component, for program purposes, is defined 
as a curved section, i.e., a section of constant curva- 
ture, a straight section, a spring hanger, a flange or 
a coupling. Although straight and curved sections, 
flanges and couplings must each be of constant outside 
diameter and wall thickness, these dimensions may vary 
from component to component. Similarly, intensive 
properties, ive., shear and elastic moduli, and density, 
may not vary within a component, but may vary from con- 


ponent to component. 


1 

See Chap. 2, 2.3 of the main body. 
See App. B, section B.94. 

See App. B, section B.95. 
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A shear area form factor of one half is built into 
the program. It may be altered at the user's discretion 
by altering one card. 

Five common end or boundary conditions, for the 
ends of the main member and the remote ends of any 
branches, may be specified by a simple coded input. 
Unusual boundary conditions may be considered by recourse 
to a more complicated input % 

Systems having rigid intermediate supports and/or 
non linear elements may not be analyzed by the use of 
this program. The number of natural frequencies found 
for a given system is a function of the significant 
figure capacity of the computer employed since the num- 
ber of significant figures required in a solution is a 
function of the characteristics of the system under 
analysis, and for any system increases with increasing 
mode number. 

Ae2 Progran Structure 

The program consists of a six section main body and 
twenty-four subroutines. If sine, cosine, sinh, and 
cosh function routines are not available in the FORTRAN 
library they must be provided by other means since, 
although required by the program, they are not explicitly 
contained therein. The function of each section of the 
main body and of each subroutine is given below. 

see App. B, section B.91. 

2See App. B, section B64, 
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Ao21 Main Body Sections and their Functions 


INPUT Control read-in of data, allocate storage 
locations for arrays and subscripted vari- 
ables. 


INVARIANTS Compute, and assign appropriate subscripts 
to, such system component characteristics 
as moment of inertia, and radius of gyra- 
iB iLoval « 


CONTROL 1 Construct transfer matrix of the system 
for the in-planelcase by calling for 
various subroutines in a sequence governed 
by the system geometry, the mathematical 
model specified, and the desired assumptions. 


CONTROL 2 CONTROL 2 performs the same function as 
CONTROL 1 but for the out-of-plane“case. 


ITERATION Apply boundary conditions to the system 
transfer matrix constructed by CONTROL 1 or2 
so as to obtain the system frequency 
determinant. Evaluate the determinant 
and utilize its sign’in controlling the 
iterative search for mode frequencies. 
Initiate output when the specified fre- 
quency range has been fully explored, the 
required number of modes found, or the 
Significant figure limit of the comouter 
reached. 


OULPUT Control output format. 
A.22 Subroutines and their Functions 


SUBSEC Subdivide components to be treated by 
lumping mass on the basis of L/D ratio 
(or ratio of radius of curvature to diam- 
éter ‘and included angle of arc, as 
applicable.) 


DISTM Construct transfer matrix of straight 
sections having distributed mass for Case 1. 


laenceforth Case i. 

Henceforth Case 2. 

See App. B, section B.94. 
4uenceforth RHO/D where RHO is the radius of curvature 
of the component. 








DISTMO 


SFIZLO 


CFIELD 


CFIELO 


POINT 


POINO 


MATMUL 


RIGID 


RIGIO 


STIFCO 


STIFOO 


HANGER 


HANGEO 


BRANCH 


Construct transfer matrix of straicht 
Sections having distributed mass for 
Case 2, 


Construct the transfer matrix of straight 
massless subsections for Case i. 


Construct the transfer matrix of straight 
massless subsections for Case 2. 


Construct the transfer matrix of curved 
massless subsections for Case i. 


Construct the transfer matrix of massless 
curved subsections for Case 2, 


Construct the transfer matrix of a point 
mass for Case 1. 


Construct the transfer matrix of a point 
mass for Case 2. 


Construct the transfer matrix of those 
components treated by lumping mass through 
the ordered matrix multiplication of the 
transfer matrices of massless subsections 
and point masses. 


Construct the transfer matrix of a rigid 
straight section having distributed mass 
for Case 1. 


Construct the transfer matrix of a rigid 
straight section having distributed mass 
for Case 2. 


Construct the transfer matrix of a massless 
rigid corner (essentially a transformation 
of coordinates at an elbow) for Case 1. 


Construct the transfer matrix of a massless 
rigid corner (essentially a transformation 
of coordinates at an elbow) for Case 2. 


Construct the transfer matrix of a spring 
hanger for Case 1. 


Construct the transfer matrix of a spring 
hanger for Case 2. 


Construct a transfer matrix from the trans- 
fer matrix of a branch and its remote end 


beh) 





boundary conditions, which permits a 
pranch to be treated as a component of 
the main member, for Case 1. 


BRaNCO Construct a transfer matrix from the 
transfer matrix of a branch and its 
remote end boundary conditions, which 
permits a branch to be treated as a com- 
ponent of the main member for Case 2. 


STAVEC Construct boundary condition vectors 
from coded input for Case 1. 


STAVEO Construct boundary condition vectors 
from coded input for Case 2. 


INVERT Invert a matrix. 


FINMAT Construct the transfer matrix of the 
system by successive premultiplication 
of the partial system transfer matrix by 
the transfer matrix of each system com- 
ponent as soon as each of the latter is 
constructed. 


FINBRA Construct the transfer matrix of a branch 
by successive premultiplication of the 
partial branch transfer matrix by the 
transfer matrix of each branch component 
as soon as each of the latter is constructed. 


DELMA Construct the frequency determinant of 
the system by applying the boundary con- 
ditions of the main member to the com- 
pleted system transfer matrix. 


eS) Program Nomenclature 
A.31 Main Body 


A Array/matrix. Transfer matrix of a 
massless subsection. 


AA Discriminant. A control in the conver- 
gence sequence. 


AAMU In INPUT; density (lbs/cu.ft.) of the 
first system component, 








AMU 


AD 


ADOL 


AFR 


AE 


AG 


AID 
AIX 


ALY 


AJT 


AMWMYK 


BB 


In INVARIANTS; average mass per unit 
length of the main nember. 


In INPUT; component density (lbs/cu. ft.). 
In remainder of programs component mass 
per unit length. 


Average outside diameter of the main 
member . 


Current computed iteration frequency 
increment. 


Average moment of inertia of the main 
member o7 


Elastic modulus of the first system 
component. 


Shear modulus of the first system com- 
ponent. 
h 


Average inside diameter of the main member; 


Component radius of gyration about the 
X axis. 


Component radius of gyration about the 
y or Z axis. 


Component moment of inertia about the 
Yor Zlavcic. 


Component polar moment of inertia. 

Total length of main member. 

Discriminant. Controls choice of mathe- 
matical model, i.é., lumped or distributed 
maSsS o 


Matrix. Transfer matrix of a point mass. 


Discriminant. Controls iteration fre- 
quency during convergence. 


lemployed in computing initial frequency increment and 
solution acceptability criterion. 

~YEmployed in computing the initial frequency increment 
and solution acceptability criterion. 
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BDO 
Bail 
Bit 
BR 


cc 


CLA 


CLY 


CLZ 


CR 


CTX 


DI 


Dr 


Dom 


Initial frequency increnent. 

Component density (lbs/cu. ft.) 

Number of last mode found. 

Initial solution acceptability criterion. 


Discriminant. Designates final iteration 
during convergence. 


Linear spring OO MENU ah of a hanger in 
the x direction. 


Linear spring cop Seane of a hanger in 
the y direction. 


Linear spring constant of a hanger in 
the z direction, 


Current value of the solution acceptability 
criterion. 


ROCHE Se constant of a hanger about 
the x axis. 


Torsional Spr ane constant of a hanger about 
the y axis. 


Torsional peeing constant of a hanger about 
the Z axis. 


Component outside diameter. 

Discriminant. Controls read-in and stor- 
age of the elastic and shear moduli and 
density of components. 

In INPUT; component wall thickness, 

in remainder of program; component inside 
diameter. 

Component wall thickness. 


Current frequency increment during con- 
vergence. 


1-6No coupling is permitted in action of hanger. For 
treatment of hangers in which coupling occurs see 


App. B, B.9. 
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E Component elastic modulus. 


BE Average elastic modulus of the main 
member. 

FFac Shear area form factor. 

FR First mode frequency of a fictitious 


straight beam having fixed ends, physical 
properties and transverse dimensions 

equal to an average of the characteristics 
of the main member to which it is equal 


in length. 

G Component shear modulus. 

GG Discriminant. Control in test for approach 
to significant figure limit of computer. 

HDEL Reference value of the frequency deter- 
minant (significant figure limit test). 

HM Number of highest mode sought. 

i Subscript. 

ite Subscript which associates a branch with 
the angle at which it joins the main 
member. 

IOP Discriminant. Controls analysis to pro- 


vide Case 1 and/or Case 2 mode frequencies 
as specified. 


JK Discriminant. Controls solution sequence 
during Case i or 2 analysis. 


KK Discriminant. Controls output to provide 
a print-out that an analysis has been 
terminated due to significant figure limi- 
tations when appropriate. 


L Iteration number. 

LM Number of current iteration, 

M Number of mode currently being sought. 
MH Number of highest mode found. 

N Subscript which associates a branch with 


its remote end boundary condition. 
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NBR 
NM 
NMBR 
NN3BR 


NS 
OMeGaA 
OMEGAC 


OMEGAG 


OMEGAT 


OMEGAM 


Number of branches in a system. 
Discriminant. Controls output. 

Total number of system boundary conditions. 
Discriminant. Assigns a component to 

the main nember or to a branch, and if 

the latter, indicates its position in 

the branch. 

Total number of system components, 

Current iteration frequency. 


Convergence reference frequency. 


Upper limit of frequency range of investi- 
gation. 


Investigation starting frequency. 
Mode frequency. 


Sum of the numerical values of the spring 
constants of a component. 


Value of the frequency determinant used 
as a reference during iteration. 


Angle at which a branch joins the main 
member. 


Reciprocal of component flexural rigidity. 
Component radius of curvature. . 


Discriminant. Controls consideration 
of rotational inertia. 


Discriminant. Controls consideration of 
Shear deflection. 


Coded boundary condition for Case 1. 
Non-coded boundary condition for Case 1. 
Coded boundary condition for Case 2. 
Non-coded boundary condition for Case 2. 


Discriminant. Controls significant figure 
CWE Or fe, 
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THuTta Included angle of are of curved com- 
ponents in degrees in input, radian 
measure in remainder of program, 


THETR Included angle of are of curved com- 
ponents in degrees. 

TL Total length of a component. 

U Array name. Component transfer matrix. 

UU Array name. <System transfer matrix. . 

V Array name. System or branch transfer 
matrix. 

VV Array name. Branch transfer matrix. 

vA Number of subsections of a component. 


4.32 Subroutines 





RHOD Absolute value of radius of curvature of 
component. 

RA Component L/D ratio. 

RD Component RHO/D ratio. 


b. DITM and DISTMO 
ARA Component cross sectional area. 


AREA Component equivalent shear area. 


co02 Array elements and/or computed constants 
C03 used in evaluating array elements. 
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3K Component flexural rigidity. 
W Component torsional stiffness (DISTMO). 


Gq ovsibelyd) eyoyel Selle: 


ARA Component cross sectional area (SFIELD only). 
ARGA Component equivalent shear area. 
bi Length of subsection. 


d. CFISLD and CFIELO 


ARA Component cross sectional area, 

ARZA Component equivalent shear area (CFIELD only). 
Cr 

SP 

Fi Array elements and/or constants used 

F3 in computing array elements. 

F5 

F6 

W Component shear stiffness (CFIELO only). 
Pat Included angle of arc of subsection. 

RD Absolute value of radius of curvature 


(CFIZLD only). 


W Component longitudinal stiffness (CFIELD 
only). 
RHOD Signed component radius of curvature 


(CFIELO only). 


SD Discriminant. Controls consideration 
of shear deflection (CFIELO only). 


V Factor. Imparts sign to certain array 
elements. 


e. POINT, POINO, PIGID, and PIGIO 

AM Total mass of subsection or component. 
ff.  obireco and SiLroo 

V Factor, Imparts sign to certain array 


elements, 
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&e BRANCH and BRANCO 


VV Branch transfer matrix (BRANCH) 
UV Branch transfer matrix (BRANCO). 
SVB Array name. Remote end boundary condi- 


tion vector (BRANCH). 


SR Array name. Remote end boundary condi- 
tion vector (BRANCO). 


Al 2 Array names. Submatrices of branch 
transfer matrix. 


R Array name. Matrix product of R2 ° ni~t 
in first usage, spring matrix of branch in 
second usage in BRANCH. 


S Array name. Spring matrix of branch 
BRANCO. 

G1, G2 Coordinate transformation matrix elements. 

U Array name. Transfer matrix which repre- 
sents a branch as a component of the main 
member . 
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APPENDIX B 


INSTRUCTIONS FOR PROGRAM USE 


B.1 General Remarks 

Appendix B is intended to provide complete instruc- 
tions for using program VIPIPE. To that end, since for 
problem programming purposes, most system components must 
be identified with respect to a coordinate reference, a 
description of the reference system employed is given first. 
It is followed by a summary of data required, a description 
of the component identification system, an explanation of 
the nature of a state vector in terms of whose elements 
the system's boundary conditions, and its displacements 
and internal force systems are described, and a description 
of the data input format. The appendix concludes with 
some comments on the options available in the progran. 

For the purpose of setting up a problem for solution, 
the longest continuous run of pipe, or the one from which 
all branches emanate if these are not one and the same, 
must be designated as the main member. A system component 
is defined as a section of pipe, straight or of constant 


curvature, an elbow, a spring hanger, a coupling? 


or a 
flange. Diameter, wall thickness, elastic and shear moduli, 
and density may vary from component to component, but may 
not vary within a component. Since both torsional and 


linear spring rates (for hangers) apoear internally as 


lsee App. B, section B.95. 
Bel 








elenents of a matrix used in obtaining matrix products, 
they must have finite values. Thus perfectly rigid 
hangers may not be treated exactly, but they can be 
approximated by using very large spring constants. 

B.2 Coordinate System 

In accordance with the procedures given below, a 
local coordinate system is constructed at each point of 
interest in the system. These points are usually the 
ends of each component, the beginning and end of the main 
member, and both ends of each branch. The coordinate 
systems are established with the piping system quiescent. 
Once established, they remain fixed in space. 

In establishing a coordinate system, one end of the 
main member, henceforth called the left end, must be chosen 
as a Starting point. 

A right-handed coordinate system is used whose x 
axis, always directed away from the left end, is the local 
tangent to the central axis of the main member. The origin 
of the coordinate system lies on the main member central 
axis, and the xy plane coincides with the quiescent plane 
of the system. The Zz axis is always directed in the same 
direction, and this direction shall be called downward, 
The direction of the local y axis is chosen so as to com-= 
plete the right-handed system. | 


entral axis 





Fig. B.2 Coordinate System 
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The starting point, for reference purposes, of a 
branch is taken as its far end, 1.e., the end farthest 
removed from the intersection of the branch with the main 
member. The plane of the branch coincides with the xy plane. 
The positive Z axis of the branch and main member are 
parallel and have the same sense. The y axis is oriented 
so as to complete a right-handed system. 
Be3 Branch Joining Angle 
The angle at which a branch joins the main member is 
measured counter-clockwise in degrees from the x axis of 
the main member as it is oriented immediately before the 
intersection, to the x axis of the branch as it is 
oriented at the point of intersection. 
B.4Y Information Required to Program Problem 
B.41 General 
&@.o number of highest mode sought. 
bo frequency range over which investigation is 
to be performed (initial and final frequencies 
in radians per second). 


c. whether in-plane out-of-plane or both in-and 
out-of-plane mode frequencies are sought. 


ds. whether employment of a distributed or of a 
lumped mass model is desired. 


ee whether it is desired to consider or ignore 
Singly or together the effects of shear 
deflection and rotational inertia. 

f. boundary or support conditions at the ends 
of the main member and at the remote end of 
each branch. 

Go number of components in the system. 

he number of branches in the system. 


B-3 








B.42 All Components sxcent Hangers 


Bo 
Do 
Co 
do 


Eo 


outside dianeter (inches). 
wall thickness (inches). 
density (lbs./cu. ft). 
shear modulus (psi). 


elastic modulus (psi). 


B.43 Straipht Sections 


Bo 


length (inches). 


Bot4+ Curved Sections 


Bo 


De 


Be45 Hangers! 


Ao 


De 


Eo 


f. 


Signed radius of curvature (negative if 
center of curvature lies on the negative 
y axis, positive if on the positive y axis). 


arc central angle (degrees). 


linear spring constant in x direction 
(lbs; 7ins). 


linear spring constant in y direction 
(Tbs ./lne) 


linear spring constant in z direction 
(lbs ./ins)s 


torsional spring constant about x axis 
(in. =-Lbs< /rad- ) C) 


torsional spring constant about y axis 
(Ane=Lbso/radede 


torsional spring constant about z axis 
(in.-Lbs./rad.). , 


Bed Component Identification System 


Components are identified by two numbers. The first, 


a component sequence number, is used to order data cards 


containing component extensive and intensive property data. 


lSee app. B., section B.96, Chap. 2, section 2.2a. 


Bah 





The second, a branch component identification number, is 
read in as data and serves to indicate whether or not a 
component is a member of a branch and if so, the position 
it occupies in the branch. 

In addition, all branches, if any, are assigned a 
sequence number which is used to order data and/or data 
cards pertaining to branch boundary conditions and joining 
angles. 

The methods of assigning these numbers are given 
below. 

B.51 Component Sequence Number 

All system components bear a component sequence 
number. The number one is assigned to the leftmost com- 
ponent of the main member. Each succeeding component 
encountered in tracing the main member is assigned a 
number one greater than that assigned to its lefthand 
neighoor. If and when a branch is encountered its com- 
ponents are also assigned sequence numbers beginning with 
the component farthest renoved from intersection of the 
branch vith the main member. That component is assigned 
a number one greater than that assigned to the last com- 
ponent of the main member encountered prior to, and abut- 
ting on, the junction. When more than one branch joins 
the main member at the same point, the user is free to 
choose the order in which desired to take up the branches 


for the purpose of assigning sequence numbers to their 
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components. Having established the order, the first 
branch is treated as if it were the only branch. The 
remaining branches are treated in exactly the same nanner 
except that the component of the preceeding branch having 
the highest sequence number, i.e., the one abutting the 
joint, assumes the role performed by the last encountered 
main member component in treating the first branch (see 


lige. Bis-5 i.) « 


Bem end 






©-Spring hanger 


Fig. B51 Component Sequence Numbers 
B.52 Branch Component Identification Number 


Assignment of a component to a branch is indicated 
by the value of a fixed point constant (NNBR) coded as follows: 
&e NNBR=1 if component is first member of a 
branch, i.e€., the component 
farthest removed from the inter- 
section. 


be. NNBR=2 if component is neither the first 
nor the last member of the branch. 


c. NNBR=3 if the component is the last 
member of the branch. 
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dad. NNBR==t if the branch has only one cone- 
ponent. 


é. NNBR=0 all main member comvonents. 
B.53 Branch Sequence Nunber 
Branches are numbered sequentially as they are 
encountered in tracing the main member. If two or more 
branches join the main member at the same point, the branch 
whose components have been assigned the lowest sequence 
numbers is assigned the lower or lowest branch sequence 


nuuber (see Fig. B.53)- 






branch 7 2 








branch # 3 


= 


Fig. B.53 Branch component identification 
and branch sequence numbers. 


B.6 Boundary Condition Codification 

A boundary condition of the main member, or of a 
branch if any, consists of six parameters which, in the 
program, constitute the elements of a column matrix called 
a state vector. Five comion boundary conditions, those of 
a fixed, free, pinned, propped, or a roller supported end 
may be indicated by a code given below which the program 
can interpret, and from which it constructs the required 


State vector. 
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2.61 Boundary Condition Code 
Codable boundary conditions are indicated by the 
below listed values of the floating point constant SV: 
a. SV=1. fixed end 
be SV=2,. free end 
c. SV=3. vinred ena! 
d. SV=4,. propped end2 
e. SV=5. roller supported ena? 

when a boundary condition is other than one of the 
five expressable by this code, the user must construct the 
necessary state vector and cause it to be read in as data. 
This situation will almost always obtain when both Case 1 
and Case 2 mode frequencies are sought for the same system. 
It is apparent, for instance, that a constraint which fits 
the definition of a prop when considering Case 1 vibrations, 
may not be describable by any of the codes available when 
Case 2 vibrations are considered. 

In constructing a state vector it is only necessary 
to know whether the value of a constituent parameter is 
zero or non-Zero. The code for a parameter having a value 
of zero is 0., while that for one having a non-Zero value 
Ars! Leys 

The constitution of the state vectors for Cases 1 


and 2 are given below. 


i Bs App. B, B.64 Catalog of Coded State Vectors 





B.62 In-Plane State Vector (Case 1) 
The vector elements in order are: 
u displacement in the x direction. 
v displacement in the y direction. 
Vp slope, in the xy plane, with respect to the x axis. 
My moment about the Zz axis. 
Vy Shear in the y direction. 
N normal force in the x direction. 
B.63 Out-of-Plane State Vector (Case 2) 
s twist about the x axis. 
T torque about the x axis. 
-w (negative of) deflection in the z direction. 
W Slope, in the xz plane, with respect to the x axis. 
My moment about the y axis. 


Vz shear in the z direction. 


Sa axis 


Fig. B.63 Vectorial representation 
of state vector parameters 
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B.64 Catalog of Coded State Vectors 


a» Case to 









tate 
Vector 
Parameter 


Vector 
Parameter 


lo 
0. 
ie 
Oo 
1. 
0. 





B.7 
Bei ¢ 1 


problem. 


Format of Data 
Problem Statement 


The first data card contains a statement of the 


The card has ten fields under the control of the 


following field specifications: 


312 , 4F3.0 , 2F8.0 , F2.0 


Field names and the information each field conveys 


are given below followed by an explanation of codes used, 


Field 


1. 
Co 
30 


Jo 


10. 


Name and Content 


NS 
NBR 
IOP 


HM 
AMYK 


SD 


RI 


OMEGAG 


OMEGAT 


GDaBC 


number of components. 
number of branches. 


code calling for in and/or out-of-plane 
solution 


number of modes sought. 


code calling for the use of a lumped or 
distributed mass model in analysis. 


code calling for inclusion or neglect of 
shear deflection effects. 


code calling for inclusion or neglect of 
rotational inertia effects. 


upper limit of frequency range of investi- 
gation (radians/sec.). 


lower limit of frequency range of investi- 
gation and initial interation frequency. 
If a zero is placed in this field the 
search for mode frequencies will begin at 
001 radians per second. 


code calling for print out of the frequency 
and value of the frequency determinant for 
each iteration, the number of subsections 
assigned each component treated by lumping 
mass, the boundary condition state vector 
for both ends of the main member and of each 
branch, the number of components in the sys- 
tem, the frequency range of investigation, 
the highest mode sought, and all of the input 
data. 
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Codes used in the first data card are as follows: 


ae LOP=0 both in-and out-of-plane mode fre- 
quencies sought. 


b, ILOP=1 in-plane frequencies sought. 

Cen LOr=Z2 out-of-plane mode frequencies sought. 
de aAMYK=1 lumped mass model desired. 

e. AMYK=O distributed mass model desired. 

f. SD=0 consider shear deflection. 

ae opal neglect shear deflection. 

he. RI=0 consider rotational inertia, 

i. RI=1 neglect rotational inertia. 

j.- GDABC=0 print solution only. 

ke GDABC=1 orint solution and input data. 


1. GDsASC=2 print solution, input data, and 
grapn data. 


B.72 Boundary Conditions 

Boundary conditions for the in-and out-of-plane case 
are read in separately. Codes outlined in B.61 apply to 
both cases. When both in-and out-of-plane mode frequencies 
are sought, the boundary conditions for the in-plane cases 
are read in first. 

The field specification for the first boundary condi- 
tion data card is: 

Cal 30 

The first N+2 fields of this card must be filled where 

N is the number of branches in the system. The coded repre- 


sentation of the main menber's left and right end boundary 
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conditions are entered in the first and last fields used 
respectively. The codification of the far end boundary 
condition for each branch is entered in the intermediate 
fields in accordance with the branch sequence number (B.53). 
When a boundary condition is not describable by one of the 
codes given in B.61, a 0. must be entered in its data field. 
For each zero entered in the first boundary condition data 
card an additional data card is required on which is entered 
the actual state vector as given in B.62 or B.63, codified 
in accordance with the last paragraph of B.61. These carast 
are placed immediately behind the first boundary condition 
data card their order dictated by the sequence number of 
the branch to which they apply (B.53). 
B.73 Branch Joining Angle 

When the system includes branches, theangles at which 
the branches join the main member are entered in one or two 
data cards as required whose field specification is: 

20 F 7.3 

For an N branch system the first N fields are used. 
The angles, in degrees, are entered in successive fields 
in accordance with the sequence number (B.53) of the branch 


to which they apply. 


B.74 Component Extensive Properties and Branch Component 
Identification 


A data card is required for each system component. 


The order of these cards is dictated by their component's 


lrielg specification 6F2.0. 
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sequence number (B.51). The card field specification is; 
3P 5.2 , MA2P%Gs2 , (GPs eee 

Field names and their data content are given, in order, 

below: 

D outside diameter (inches )4 

DI wall thickness memes = 

RHO radius of curvature with sien (B.44) 
(inches), Note: Applies to curved com- 
ponents only. A zero (0.) is entered in 
all other cases. 

AO length (inches). Notes applies to 
straight components only. A zero (0.) is 
entered for all other cases, 

THSTA are central angle (degrees). Note: Curved 
components only. A zero (Q) is entered in 
all other cases. 

CLA 

CLY 

CLZ \spring constants (Bo4sy 

CTX 

CaN. 

CrZ 

NNBR branch component identification number (B.42). 

B.75 Intensive Properties of Components 
The first intensive properties data card has four 


fields whose specifications are 


feo, Cheol career 2.0 


i-20nese fields are left blank in cards pertaining to 
spring hangers. 
These fields are left blank except in cards per-= 
taining to spring hangers. 


Bei4 





The field names and their data content are given below. 


AAMU density of first component (lbs./cu. ft.). 


AE elastic modulus of first component (psi). 
AG shear modulus of first component (psi). 
DD code. see * below. 
% When the intensive properties of all components (spring 


hangers excluded) are identical, the DD field is left blank 
and the first intensive properties data card becomes the 
only intensive properties data card. If DD is non-Zero, 
the data on the first card is interpreted as pertaining to 
the first component only and an intensive properties data 
card must be prepared for each of the remaining components, 
i1.e€., the second through the last, spring hangers included. 
These additional cards are identical to the first card in 
format and field content except that the DD field is deleted. 
A blank card must be inserted for each hanger. Intensive 
properties data cards are ordered in accordance with their 
component's sequence number (B.51). 
B.8 Data Deck Assembly 

The data deck is ordered as follows: 

1. problem statement card. 

2. boundary condition data card(s). 


3. branch joining angle data card(s). (deleted when 
system has no branches). 


4. component extensive properties data card(s). 


5. component intensive properties data card(s). 
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The user may obtain solutions for any number of prob- 
lems desired by “stacking” assenbled problem data card 
decks at the end of the program deck. Note that a blank 
card must be inserted between the last card of the prosram 
deck and the problem statement card of the first problem. 
A blank card must also be placed behind the last intensive 
properties data card of the last problem. DO wOT insert 
blank cards between problems. The problem statenent card 
of the N+tfirst problem must be immediately behind the last 
intensive properties data card of the nth problen. 

BeJ9 Flexibilities of the Program VIPIPE 
B.91 Shear Area Form Factor 

A shear area form factor of one half is built into 
the program. [g] It may be altered by changing card 000580. 
B.92 End Conditions 

A properly constituted boundary condition state vector 
must have three zero and three non-zero elements. fFre- 
quently, however, an actual end condition may not meet 
this specification, e.g., an end of the main member or of 
a branch may be supported, or connected to a piece of 
machinery in such a manner as to approximate being held 


in a rubber mounting block (Fig. B.92-1). To program this 


fictitious free end—— 
a 2 
system > SS 


SENT POS 


FPi-. 92-1 Fictitious end condition 
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end condition, the user should make a good approximation 

of the elastic properties of the mounting, support, or 
connection in terms of the linear and torsional spring 
rates of the program (B.45, B.74)* then attach a fictitious, 
short (perhaps .01 in.) massless straight section to the 
end under consideration, and specify in the boundary con- 
dition for that end that the end is free. 

Through the use of fictions of (this sort, preperly 
constituted end condition state vectors may be constructed 
for any real situation with little or no loss of accuracy. 
Bo 93 Subsectioning 

Subsectioning for the lumped mass treatment of com= 
ponents is carried on in subroutine SUBSEC? , The subroutine 
is set up to divide a component into a maximum of ten sub- 
sections, The writer, through the analyses of several 
systems, found that this figure gave good accuracy of 
solution for a reasonable expenditure of machine time. 

For all components having a length-to-diameter ratio 
greater than six, the number of subsections is computed 
as two times the value of a parameter HM which appears 
in the subroutine. The user may control its value by 
altering card 004720 in the subroutine. 


see App. A, sections A.22, and A.32. 
1.€., treat the mount as a hanger. 





4 





B.94 Values of the Frequency Determinant and Solution 
acceptability Criterion 


The freauencies for which the value of the frequency 
determinant is zero are the natural frequencies of the 
system which it describes. In light of the number of sig- 
nificant figures available in a computer and of the effects 
of round-off, however, insistance that iteration and con- 
vergence operations continue until the value of the 
frenuency determinant actually became zero viould require 
an excessive expenditure of ma@emime Cine., Instead, itera— 
tion proceeds until the sign of the value of the frequency 
determinant changes, then a convergence routine is carried 
out until the difference between the values of frequency, 
for which the corresponding valve of the frequency determinant 
have opposite signs, is less than a svecified value, CR. 

The result of a linear interpolation between those fre- 
quencies is accepted as a natural freaquency of the system. 
The specified frequency difference CR is termed the solution 
acceptability criterion. 

To compute a starting frequency increment, the prograa 
constructs a synthetic straight pipe, equal in length to 
the length of the main member, and having diameter, wall 
thickness, and intensive proverties equal in magnitude to 
a weighted average of those properties of the composite main 
member. The fundamental frequency of the synthetic pipe 
1s computed, and one half of its value is taken as the 
starting frequency imenenenee The initial solution 


1 
This increment may be cae by altering card 001090. 
Bel 





acceptability criterion is taken as one-twenty-fifth of 
the starting frequency increment. The initial solution 
acceptability criterion may be chenged by altering card 001110. 
The starting frequency increment and the initial solu- 
tion acceptability criterion are used until the second 
natural frequency is found, after which the frequency 
increment is recomputed, as it is after each succeeding 
natural frequency is determined, and given a value equal 
to one half of the difference between the last two fre- 
quencies found. The new solution acceptability criterion 
is then recomputed as one-fiftieth of each new frequency 
“increment. The new solution acceptability criterion may 
be changed by altering card 002800, 
B.95 Flanges and Couvlings 
Pipe couplings and flanges represent real "lumps" of 
mass mounted ina piping system. They may be treated as 
individual system components. If this is desired, a com- 
ponent sequence number and a branch component identification 
number should be assigned in the usual manner (B.51, and 
B.52), and appropriate extensive and intensive properties 
data cards (B.74, and B.75) should be included in the data 
deck. These components are treated in the program, as 
rigid sections having distributed mass, 
B.96 Hangers exhibiting Coupling 
Hangers such as that shown in Fig. B.96=-1 can be 


treated with the tools available in VIPIPE by considering 
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the mounting arm as a branch having a hanger 


Fig. B.96-1 Cantilevered hanger 
and a free end. Recourse to similar methods will permit 
treataent of almost any hanger configuration. 


B.97 systems Having Unusual Flexibilit 





When the main member of a system is essentially 
cantilevered the computed starting frequency increment may 
be so large that the first two system natural frequencies 
are not found. For analysis of such systems it is reconm- 
mended that the constant 4.73004019 in card 001080 be 
changed to 1.87510067. Similarly, if it is suspected that 
the fundamental frequency of the system is less than 
-01 rad./sec. a suitably lower starting frequency should 


be programed (B.71). 
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APPENDIX C 


FLOw DIAGRAMS 


C.1 General Remarks 

Flow diagrams are included for all parts of the main 
body, and for several subroutines. Those subroutines for 
which no flow charts are included, with the exception of 
INVZRT, incorporate, at most, two decisions in their 
Structure. As applicable they sense on SD and/or RI ( A.31) 
to change the value of one or two elements of the transfer 
matrix constructed by the particular subroutine. 

Subroutine INVERT performs a matrix inversion by a 
pivotal technique. Sufficient comment statements are 
included in the subroutine to define its structure. If the 
user's FORTRAN library includes a satisfactory single 
precision matrix inversion subroutine, INVERT may be deleted 
after suitably modifying the dimension (cards 010070, 

010080 BRANCH, and 011080, 011090 BRANCO) and call (card 
010230 BRANCH, and 011660 BRANCO) statements of subroutines 
BRANCH and BRANCO. 


C.2 Index of Flow Diagrams 


C.21 Main Body 
Flow Diagram Page 


INPUT Gas 
INVARIANTS C=5 
CONTROL IN-PLANE @a6 
CONTROL OUT-OF-PLANE C=6 
ITERATION C-7 
OUTPUT C-9 


lsee 6.3 Cot 





C.22 Subroutines 


Flow Diagram Page 
SUBSEC €=10 
MATMUL c- 2 
BRANCH cain 
BRANCO Cate 
STAVEC Cai 
STAVEO Gai 
FINMAT Ca lity 
FINBRA Catt 
DELMA Caik 

C.3. Flow Diagrams! C=-3 through ¢.14 


Note: Since CONTROL IN-PLANE and CONTROL OUT-OF-PLANz 
are identical in structure, only the former flow diagram is 
included. They differ only in the subroutines called and the 
State vectors used (4.22). Subroutines SUBSiC, MATMUL, ANVERT, 
FINMAT, FINBRA, and DzLMA are common to both cases. 


14 complete list of symbols is given in Appendix A, 
section A.3. 





» 





\ READ; Dimension 
statement 


READ: Problem 
statement 


CO ore $9 











READ: Coded Case 1 
boundary conditions 


READ: Coded Case 
boundary conditions 


J=1 


>0 









READ; Coded state 
vector SVOP(I,J) 


READ: Coded state 
vector SVIP(I,J) 


I=1,6 = 
13> a J=J41 
=0 
>0 fom 
NBR >=2 


READ; Branch joining 
ngle PHI(I) 
I=1,NBR 


£ 


READ; D(I),DI(1),RHO(I), TL(1), THETA(I) 
CLX(I),CLY(I),CL24(I),CTX(I),cTy(T), 
CTZ(1), NNBR(I} 

I=1,NS 


READ: AAMU,AE, 


AG,DD 
> 
ee READ: AMU(I),B(1),G(1) 
L= 2,NS 
=0 = ees: 
aay 
[ONESAT=,01] 
>0, Peers 





MEM ON. WIS EO ate etohein re) 


AMYK = AMYKel. 
OMETAO = OMECAT 
FFAC=.6 

L=1 

Nea 

AA= PR=CC=GC=EE=0. 
AAMT =AL=AD=AID =O. 
KK = 0 

NNN =NFR+1 

















VHETR(I) = THRTA(I) © 

DHPC GT )) = Topets IC)! 

RAMU(T) = AMU(T) 

DI(I)=DI(I) = 2.0#DI(T) 

T= 1,NS : 


GO TO INVARIANTS 





INVARTANTS FLOW DITA GKAM 


[FROM INPUT] 





50 CTRL) eCLY (CT) 40E2 (TP 
———— ~~ Crm (TD)ecivan) eer) 














THETA(1) = THETA(I)#¥41 /180. 
TL(I) = THETA (T)#ARSF(RHO(T)) 


s((F Jovd-07( I )sterd )f32.- 








AT(1)#E(T)) 
RTF ((D(T )#2 es ) 
[elavaccouei wi 


Bek. 


3 





a 









_TNVAL TANS JON. BE AN 






/HHM == HM 
iCALL SURSEC 


ee 


PR (TT) 


#0 


+ 


AL= AL+TL(T) 
AD =AD+D(1)#°L(T) 

+B =ER+E(I )e7L(1) 

AID =AID+DI(I)#TL(I) 
AAMU = AAMU+AMU( 1) L(1) 


sae 





AD= AD/AL 
AID = AID/AL 
AFR= 7/644 (AD#+4-AID#S4 ) 


FR= 4, 734#2#SQR TF ( EEH#AFR/ (AAMUSAL*4 ) 


ADOM = FR/2.0 
DOM = ADOM 

CR = FR/25, 
BDOM =ADOM 
BR= CR 





PHI(II) =PHI(IT)*T/180. 


II= 1,NBR 


<0 70 


[Mn = 1 2G 


TO CONT ROL_IN PLANE 


TO CONTROL OUT OF PLANE 








Creer os! PLA'E FT OW DIEGPAM 


[Fro (Avie 
i 
St es Spee af ’ a eth peso SS eg 
[svz(3,1) = SVIP(J,2) ~ -SL<auap t=. - >{CALL STAVEC| 
{ =e Se Soo ee 
ee J | 


i , 
SVE(J,1) = SVIP(J,NMBE ) }=——< oe Su, avin | 


sin) BS 
eee | : - ea 


— 









=2G 1p >O 


ALL RIGID CALL DISTM CALL STIFCO (ry OlCALL CFIELD 
npr (1 p22 CALL FINBRA 


[1 TrERATION] 





ITEKATION FLOW DT ATHAM. 
A [CALL DELMA | 


{determinant contructed by DELMA~ 


IDEL= DEL(L) 


GG=2. [Numerical value ofmcae fagapencs| 
B 
OMEGA (L+41) = OMESA(L)+2.#CR 


>0 sh 
<< 
is 





DELL) =-sle= 1. 
%O 








>O 


ns >toprer=DEiC 
=0 
PDEL = DEL( T= 









MEGAM(M) = OMEGA (L)+( OMEGA (Le1 )-OMEGA(L) )#ABSF(DLL(L 


ABSF(DEL(L) )+A® SF(DEL(L=1 


DELC = DEL(L) 20r7. >0 < 
OMEGAC = OMESA(L) a> | BR = 0. | 
=O 
<O 


OMEGAC = OMEGA (L=1) 
=0 AA= 1 e 


(A)}<—fOMEGA(T+1) = OMEGA (L )+DOM <iraca(1)- ovEca SA We Wel] | NAR BB=1 
) Tourn 50 DELC = DEL(L=1 


JB. x 22KPB 
TA aa 
[apex 3.] ‘ 








Pan 


OMEGA (I+1) = OMEGA (L)=( OMEGA(L)-OMEGAC ) 





>0 







ae OMEGAC +( OMEGA (L)-OMEGAC )#APSF(DELC <0 oe 
AESE(DELC ICAREP DET TT OMECGA(L)=-OMEGA(L 1)-CR 


bac~o +4 |? 
AA= BB =CC =GG=0 
DOM = ADOM 


ADOM = (OMEGAM(M-1 )-OMEGAM(M-2) )#.8 
CR= ADOM/50. 


OMEGACL+4)=OMEGAM(M-1) + 2.%C 























OMEGAO = OMEGA(I+1 | 


—$— 


. ae OMEGA (L+1) = OMEGA (T, =. 9#D0r" 
GOTO Aldvess loo | DOM = . 5#DOM 


CON IROL. OUT -OF > PLANE —~ 





Cart 








SUBROUTINE MATMUL FLOW DIAGRAM _ 

















C(J,K)=C(J,K)}+B(J,L)#A(L,K) 
L=1,6 K=1,6 J=1,6 







WUT = cle) 
K=1,60 J=1,6 





pane Ky = Bia KO(T, mee 
iy K=1,6 90 


RETURN] 


c-8 





OUTPUT FLOW Din . AM 


dewey ae 





\BrINfifeadingy- = 2 


RINT: In-Plune, 
frequencies 





PHINT?: 








PRINT: Ma Ke. 
model used 


PRINT; Input ssi ani 


OMEGAM(M) =O. 


M=1,MH 


OMEGA(L)=0. 
L=1,LM 


DELC = OMEGAC = 





) Add ess 450 of 
INPUT loli Soon 


Caw 


20 







ADOM = BDOM 

CR = BR 

OMEGA (1) = OMEGAI 
= OMEGAI 


GO TO Address 998 of 
CONTROL OUT-OF-PLANE 


Out-of-Plan 


ay, 


requencies -— 


PRINT; Subsection data 
raph data,input data 
ee 


MN = AA= BB=CC =GG=KK =0. 








£2 0} — rine] 


Z=4e RETURN 








RUE TAIN 1G PAWEL Ia ives 8 PM 


j=] ee 
fv = yj eta) = WET KNOT Tes] OT 1) 
I {141,58 4 =1,3 J=1,3 
Pea) . 
i. aie tae ocn 
n TMIKJ=O. - 
! “if SO t=1,6 aoe 
| pee Eas) | 
obi .. | ee esa : 
as TCI, 1) =ait2, 8) = sea yet. 
_ i... AS A= (ob ac 5, 6) sale 
[M= M+] | 
Wa ;e)= 80,1) 
| W(B,L) = R(2,2) 
| 45,1) =k(3,L) 
— er eee! 
CRETIRT | 





CALL INVERT 
(invert R1) 









R(J,K)= R(J,K)+R2(J3,L)#R1(L,K) 
L=1,5 K=1,5 J=1,5 


SP= SINF(PHI) 
CP = COSF( PHT ) 


Construct third order coordinate 
transformation matrices, Sl] ard 92 





R1(7,K)= R1(7,K)+92(7,7,)4R(L,K) 
Bel,5 Ke1,5 FH1,5_ 


R(1,K) = 0. 
T=1,5 K=1,5 





Catt 








S"FROUTINE PRARCO FLOW | DIAGKEM 


wed | a a 


Reorder bourdury condition state E ) 
vector such that the first three J 
elemerts pertain to displacements, 

the last three to forces __ 





"S(J,K)=S(J,KHR1(J,L)#G1(L, 1K) 

hs ae, a ee. ee . ee a Keil, $ J=#1,5__ 

[seorder branch transfer matrix 

lag make it compatible with the { 
reordered state, vector. : sae ; 














3 [ u(1,1)= 0(2,2)=0(3,3) 21. 
. u(4,4)= U(5, 5)=u(6,6)=1. 
=O SR (Nn MT) U(2,l)=S(1,1)  u(5,4)= 8(2,3) 
ae ; U(2,3)=8(1,2) U(6,1)=8(3,1) 
i U(2,4)=8(1,.3)  U(6,3)= 8(3.2) 
U(5,1)={S(2,1) U(6,4)=S(Sy3) | 
U(5,3)= S(2,2) 


me 


"RETURN ] 


R1(1-3,K)= UV(1=3,K) 
R2(1-3,K) = UV(4-6,K) 
a ds Sains, 







CALL INVERT 
(Invert Ri 






ue) =O. 
J=1,5 K=1,5 
R1(J,K)= R1(J, eG. L)#R(L,K) 
| ‘L=1,3 K=1,3. 5=1,3 | 


Cale 








SURFONTINE STANEC FLOW Dik AN | 









= [ROTET I= HCC 5,8) = BECO N= 2 


ROE NI= FOS, N YS BFS TST | ——- > Ft] 


PC(L,V)= 0. 
L=1,6 






BO(2,N)= BC(5,N)= BC(6,N)= 1. 


BC(1,N )= BS(3,N)= BC(4,N)= Le 


RC(2,N)= BC(4,N)= BC(6,N)= 1. PETIRN 


= efRO(AgN > BO(E Y= BLO, fee 


Ca 








_SUPRQUTINE FisWAT FLOW | DIAGKAM__ 


jee! lf VUT,K) = UCT, K) +} ~- —- > Eerurw] 
ce (i= 1,6 K= 1.6) 


TWIT, K)= 00 p= WC) = VIS, K)FU(T, L)#UU( LK) b — ——[ Rerun] 
wel, x=1,6 | = te ? k= 9 J=1,607 | 


[ V(a,«) = U(a,K) [-—-——-LEeninN] 
Rac aie pe Saas 1,6 





wa K)=0. --—--- me: J=V(T,K 4U( 0, L)4VV(L,K) + RETURN 
J=1,60 K=l, 6 apne 6 K=1,6 = es 

















ee On sly 


~ 
x 


aes ‘ Ss Sar - 








APPENDIX D 


PROGRAM LISTING 


D.1 Index 


D.11 Main Body 


Section Sequence Numbers Page 
INPUT 000040-000760 D-2 
INVARIANTS 000790-001170 D-3 
CONTROL IN-PLANE 001200-001810 D=4 
CONTROL OUT-OF-PLANE  001840-002450 D=5 
ITERATION 002470-003260 D-6 
OUTPUT 003280-004660 D=-8 
D.12 SUBROUTINES 
SUBSEC 004690-005140 D=10 
DISTM 005170-005710 D=-11 
DISTMO 005740-006280 D-12 
SFIELD 006310-006640 D=13 
SFI“LO 006670=006910 D-14 
CFIELD 006940-007370 D=14 
CFIELO 007400-007890 D=15 
POINT 007920-008090 D=-16 
POINO 008120-008310 D-17 
MATMUL 008340=008900 D=17 
RIGID 008670=008649 D-18 
RIDIO 008930-009160 D-18 
STIFCO 009190-009410 D=-18 
STIFOO 009440-009660 D-19 
HANGER 009700-009850 D=-19 
HANGEO 009880=-010030 D=20 
BRANCH 010060-010730 D-20 
STAVEC 010760-011040 D=-21 
BRANCO 011070-012210 D=22 
STAVEO 012240-012480 D-24 
INVERT 012510-013220 D=25 
FINMAT 013250=013380 D=26 
FINBRA 01341 0-013540 D-26 
DELMA 013570-013750 ° D-27 
De Listing D=2 through D-27 








Anan 


PRCURAM VIPT PRC 


INPUT 


NM=0 
JJ=i 


HSOODIMENSIUN AJT(SC), AJ(SO)2C2C)sATX(SOISATY(SOPSPHI(2CIsSVI (O91), 


1500 


1501 
1502 


1503 
1504 
1505 
1506 
1507 


1508 
15¢c9 
1510 
1511 
1512 
1513 
1514 
1515 


1516 


1517 


1518 


1519 


1520 
1521 


TSVB(6),0(50) C1 (SC) ,TLESCI,2(50),G(501,E(50), 


AMU(50), THE TA(50),RHO 


2°50) ,CMEGAM(20),CEL(SCO) pCMECAECC)s A(6,b),P (65h) ,Ul 096) sCLX(5C)5C 
BLY(5C) CLZG50),CTXCSO) CTYC EO CIZISCI pVEb Ede VV( 6,6) LUCE rb) SVE ( 
WGel)eSV(Z2), SVEP CEs 220 ySVOC C219 SVOPL 6922) +NNBRE SO) sBAMUL-SC) » THETR ( 


520) ,0T(5C),P(20) 


REAL 15CCeNSsNEReTIOP,HM, SMYK SD, RT pIMECAGSOMEGAT ,GDABC 


FCRMAT(31Z,4FZ.C,2FE.0eF 2.0) 
NYBR=NBR4+2 

IF CICP-1)1561,15601,15C6 
REAL 1502,(SV(1I),1=1,NMBR) 
FOCRMAT(22F3.0) 

cC 1505 JjJ=),NMBR 
IF(SV(S1)}1503,1°02,1505 
READ 1504,¢SVIP(I,J),1=1,6) 
FOCRMAT(CFZ.0) 

CCNTINUY 

IF (10P)15C6,150€,1510 

READ 1507,(SVC(1),T=1,NMER) 
FCRMAT(22F3.0) 

CC 1910 J=l,NKMBR 

TF (SVC(J))1508,1508,151C 
REAL 1509 (SVCP(I,J),1l=156) 
FCRMAT(GF2Z.0) 

CONTINUE 
IF(NER)1572,7514,1511 

RFAD 1912, (PHI({),1=1,N@R) 
FORMAT (1CF7.3) 

OREAD 1514,(C(1),CI(1),RHC(I), 
TIXCL)sCTY(1) »CTZ¢1) sNNPR (I) oI 
FORMAT (3F5e2y2FE.2,0F7.2,12) 

READ 1515 , AAMU,AE,AG,CL 
FORMAT (F7.3y2E10.2,F2.C) 

IF(00)1516,1516,15178 
OG 1517 1=t,NS 
AMUC1T)=AAKU 

EC LT)=AE 

GUE) =AG 

GC iC 1520 
AMU(1)=AAMU 

ECT)=AE 

G(1)=AG 

READ 1579,0AMU(I),ECI),GUL),1=294S) 
FORMAT (F7.322E10.2) 

LF (CMEGAT)1520,1570,1521 
GMEGAI=.01 

CMEGA(1)=CMEGAL 


TL(1), 
=1,A\S) 


D-2 


THETACE) sCLX(TIoCLV(T ),CLZC011),€ 


OC0OCz0 
OCcCO30 
ooccKyc 
occcse 
ccocéo 
occc?70 
OoOCcAC 
ococgc 


-000100 


acol1ic 
00012C 
ocol3c 
000140 
oco1sc 
Gooléc 
ocol7c 
0co180 
000190 
o000c0c 
000210 
000220 
000230 
DO0zuC 
acozs50 
000260 
000270 
occ280 


000290 


000300 
000310 
000320 
000330 
0C0340 
000350 
000360 
0c0370 
0003e8C 
00039C 
o0cuce 
000410 
000420 
000430 
O0Ok4O 
o00ck50 
000460 
oo0okK70 
ocouso 
OCOLS0 
000500 
ocosic 
000520 
000530 
000540 
00055¢ 





Aaan 


Were 


Jol 
102 
103 


104 


105 


106 
107 


AMYK=AMYK-1.0 
CREGAC=CMEGATI 
FFAC=.* 

L=1 

v=] 

Ah=C.0 

BR=0.0 

cc=0.0 

GG=U.0 

KK=C 

AAML=6.0 
EE=0.0 

AL=0.0 

AD=0.C€ 

AID=C.0 
NNN=NBR+ 1 

DO 1922 [=1,NS 
DV(1y=01CT) 
BAMU(T)=AMU(T) 
THETR( IT) =THETACT) 
DI(iy=0U01)-2.0*C1(1) 


INVARTANTS 
OC 10@ L=1,NS 
IF(CLXCLIV4CLY CL IV+CLZ¢ 1 P+ CTXCL)+CTVEL +CTZ011) 108,101, 1c8 
AMUCT)=AMUCT) 23.14 I592654H (CCL) *#e2-DI(L)#e2)/(6912.852.1 7812.) 
IF(RHG(1))102,103, 102 
THELAC LT) =THETAC I) #3.141552654/180. 
TLOI)=THETA(I)*#ABSF(RHC(T)) 
AJT(T)=3. 1415926544 (DIL a#4-CI(T)eeu )/2c. 
AJ(TVHAST(1I/2. 
RCTI=1T-O/CAJS(TV#E 01) ) 
AIX(TP=SCRTIFC (CCL) ##24+DI (1) ee2)75.) 
AITY(LTV=AIX(CTI/1.41b2Te 
IF(AMYK) 104,105,104 
Z(1)=1.0 
IF(RHO(T))105,10€,105 
HHM=HM P 
CALL SUBSEC (RHCOT), THETAC LI, DCL), C1),HM,Z(1)) 
HM=HHM 
IFOENNPR(1T))108,107,108 
AL=AL+TL(1) 
AD=AC+O(1T)*TLtl) 
EE=EE+E(IT)*TL(I) 


. AID=AIC+ACIC I) #Tl(1) 


108 


AAMUL=AAMU4AMU(T) #TLII) 

CONTINUE 

AD=AC/AL 

AID=AIC/ AL 

AFR=3. 14 1S9CESUs OU H( ACH RU-SID HHL) 
FREU.73Z00L019e #2 eSCRIF( EER AFR/(AAMUs AL ##4)) 
ADCM=FR/2.0 


De-5 


oceséo 


0cos70 


ocassc 
oco530 
occecc 
O006IC 
occ6/70 
0606;:0 
ocodé4ua 
00065C 
accéé0 
000670 
ocoéséa 
o0c0690 
oco7c0o 
0c0710 
oco7?2c 
0c073C 
000740 
oco7sc 
o007é6¢ 
ocda770 
000780 
0c0790 
ocoecc 
ocoe8lo 
00c82C 
000830 
000840 
0c0850 
000860 
000e70 
000820 
000890 
000900 
ocog9l10 
000920 
000930 
060940 
600950 
000960 
000970 
000980 
000990 
oo1cacd 
001010 
001020 
001030 
001C40 
ocicsc 
001060 
oo1o7c 
001c80 
001090 





Ana 


DCM=AUCI 
CQ=FR/25.C 
BCOM=ACCY 
BR=ULR 
co 109 Ll=1l,8fR 
109 PHICTILISPHIC TI) *3.141592€54/180.6 
TFCIGP-1)1%10,999,998 
110 MN=] 


CONTROL IN PLANE 


999 IF(SV(1))1011,1C21, 1011 
IG11 CALL STAVEC (5V(1),1,SVI) 
GO TO 1081 
1021 GC 1631 J=1,¢ 
1031 SVI(J,1)=SVIP (J, 1) 
1641 TECSV(NMBRI)1C51,1061,1052 
1651 CALL STAVEC(SV(NMBR) 2 1sSV0) 
GO TL 108) 
1061 Ul 1071 J=1,6 
1071 SVE(J11)=SVEP (J, NMER} 
1081 JK=C 
IF(NPR11000,100C,1085 
1085 (CC 1087 N=2,NAN , 
IF (SVEN) }1087,1087, 10R6 
1086 CALL STAVEC(SVIN)NySVIP) 
1087 CONIINUE 
1000 N=2 
f=) 
00 10C1 I=1,NS 
PCTV=CLACTVFCLYC1I+CT2 (1) 
TF(P£1)137) 38,37 
37 CALL HANGER(CLA(I),CLY(1),C1Z019,U) 
GC TC 62 
38 TFCAMYK)39,41,39 
39 LFCREOC(T)1N0,42,4C 
WO TF(2(1)943,45 743 
41 LFERHC(T) 140,48, aC 
U2 TFC ZC1)24E, U7 UE 
W¥3 CALL CFIELD (RO1),Z(1)sG( 1), SD RPGC)» THETA(T) DUE), DICL) ELD, A) 
Wh CALL POINT (AMU(T),ALYET) ¢CMEGA(L),TL(T) p20 1) ,R198) 
CALL MATMUL (4,9,Z(0)9U) 
Counc 62 
WS CALL STIFCO (TRETACT)»RHC(1),U) 
GO TO 62 
HOOCALL DISTM (AMUCI),ATV(L) -GMEGA(L) sTLOT) RUT), 001) sDIC 1) GC 1),S09R 
VI,ECL),FFAC,U) 
GO IC 62 
WT CALL RIGIC CAMUCT), TLOT) sCMEGA(L),ATYCI) QT eL) 
GO TL 62 
U8 TFCZC1))89,47,49 
WO CALL SFIELD (CUCL) eTLCLD RELI eZ) 9301), S0eD( 0) ECT)» FFACZA) 
GO TC un 


oclice 


001110 
001120 
0o115¢ 
OCTI4O 
001150 
001160 
oolic 
oclisc 
0C1190 
ocizc6c 
o00121¢ 
001226 
0C1220 
oC124¢ 
oo1se 
001260 
001270 
0c128C 
001z29¢ 
0C1300 
001310 
001420 
00133¢ 
OCt2uC 
0C1350 
001360 


001376 


001380 
001390 
0c1400 
001410 
001420 
00171430 
OOV4LC 
001450 
001460 
001470 
OO14EC 
001490 
o0c1500 
001510 
004520 
001530 
0C1540 
0c1ss0 
001560 
001570 
001580 
0C1590 
001600 
oclé10 
001626 
001630 








OAD 


62 
63 


64 


65 


66 


67 
68 


1001 


UF wh 


SOCALL CFIELQCAJT( LD, RU) +200) 2G01),S0eRHG(1), THETA( I) sOCE ),OL(1). FF 
TAC, A) 
7 CALL PCINC CAMUCT),AIX(1),ALYO€1),CMEGACL),TLCI),Z01) eRI,B) 


TFONNEROL9965 763565 

CALL FINMAT (L,ULU,I.V) 

UGNOy J=1,6 

OC 64 K=1456 

UUl de KI=VOS,K) 

GO Ic 10C1 : 

CALL FINGRA (Ly VVyNNDR( L,Y) 
CO 66 J=106 

DON 66 K=t96 

VV Ode KIEV (J 9K) 

TFOKMRBR(1))68 167,67? 
TFCNNUR(1)-511001,68, 68 

CALL GRANCH(SVIP,VV,PHI( II) NU?) 
N=N¢] 

Il=11+1 

GC TC 65 

CONTINUE 

Go Tc 20006 


CONTROL CUT CF PLANE 


LF (SVC(1))1012,1023,1012 
CALL STAVEQ(SVC(1),1,SV1) 
@@ Te) VO "4 

NO 1033 J=leé 
SV1L(J,1)=SVCP (J, 1) 

TFC SVQ(NMBR)) 1053,1063,1C52 
CALL STAVEO(SVO(AMBRI,1,SVE) 
GC TL 1083 

DO 1073 J=1,6 

SVE (J, 1) =SVCP(J,NMBR) 

JK=1 

LF (NER) 1002, 1002, 1093 

DO 1095 N=2,NAN 

LF OSVO(N))1095,1095,1094 
CALL STAVEC(SVCIN) sNeSVCF) 
CONTINUE 

N=2 

(isi 

O00 1003 I=1,NS 
PCLI=CTX(LIFCLZ(1I+CTY(1) 
IF (PC1))15291 . 


CALL HANGEQ(CTX(T),CLZ(1),CTY(1),U) 


Gc TC 4 

IFCAMYK)3511,3 

IFEREC(1))4,8,4 

TF(Z01))62596 

CALL STIFCO (THETACTI),RHO(T),U) 
GC 1C 14 


D=5 


0C 1640 
001650 
00164C 
001670 
001680 
001690 
O017CG 
001710 
0C1720 
001730 
001740 
001750 
001760 
00177C 
90178CG 
0C1790 
0c18co 
O01E10 
001820 
001830 
001840 
00185¢ 
001360 
0Cc1870 
oo1eec 
001890 


001900 


CC1910 
0C1920 
001930 
001940 
001950 
0019¢0 
001970 
001980 
001990 
002000 
002010 
002620 
002020 
002040 
002050 
062060 
002070 
002080 
002090 
002100 
002110 
002120 
002130 
002140 
062150 
002160 





C 


Gere Tt 
8 TF (2(1})16,9,10 
9 CALL RICIC (AMUCED, TLUL) ,ALXC IL) ,EMEGACLI,ALTYOCL),AST(C ED oR 1,U) 
GC Tc mu 
TOOCALL DISTROCAMUCTI,ATK( IT), ALYET) ,OMEGACL I, TLIT) ,ASTEED eR (1) LUCIE 
VECEV GUT) sSCyRTy FRAC SI") 
GC [fc ju 
11 TFORFC(E)I4,12,% 
12 TFC2(1))13,9,123 
13 CALL SFIELO (CICTI, TLE) ,ASTOCT)D, 301) ,201),G601),S0,011),F FAC,A) 
GC IC 7 : 
TH TFECNNBbRO1)017,15,17 
15 CALL FINMAT (U,UL,I,V) 
CO 16 J=1,6 
DC 16 K=3,6 
16 UU(dsKI=zAVEOU,K) 
GO ic 1603 
17 CALL FINBRA (LIVV~NNER(II,V) 
NO 18 J=1,6 
oe 18 K=1,6 
18 VV(JsKI=EV(IeK) 
TFONNBRCT)) 82,159,199 
19 TFONNER(E)-2)10C3,82,°%2 
B82 CALL Bee OE SNCRs VV ePHIC IL) y Nyt) 
N=N41 ‘ 
TlslI+1 
GC Ic 15 
1003 CONTINUE 


ITcCRATION 


2000 CALL CELMA(UU,SVI,SVE,V) 
ODELC LIEV (Le TI #V( 2,2) VI 3,Z7I-VE ly SDA (C22) eV 5,1 VO Lec de Vv (2,3) eV (3 
TeDAaVET 2d eV Se 12803, 2) eV 1, 2 eV 3, 208 (2,1 RV Oe TeV (3 2) eV 2, 5) 
GO TC 61 
23 POEL=DEL(L) 
CO TC 25 
22 1F(AA)23,23,24 
23 POEL=DEL(L-1) 
GO TO 25 
24 PODEL=CELC 
25 IFC PCEL*DEL(L))26,26,27 
26 1F(AA)28,28,57 
27 BB=0 
IF (AA) 35,35,52 
28 AA=1.0 
BB=1.0 
DELC=LEL(L-1) 
CMEGAC=CMEGA(L-1) 
Go TC 59 
2900MEGAM(M)=ONECAC+ (OMEGA(LI-CMEGAC) #ABSF( CELE )/(ABSFIDELC )+ 
TABSF(CEL(L))) 
GC TC 32 
30 CMEGA(L+1)=CMEGA(L)-.980CK 


De-6 


eA 1 Sie 
Cale 
N022.CG6 
0CG2210 
CC22z20 
0C2.23C 
Oc2cuC 
Oi 2256 
0C2260 
O02z27C 
NC223C 
0c2-290 
0022C0 
NEA AWE 
90232C 
0¢2430 
0c2240 
oc235c 
002360 
0C2370 
0C2380 
002390 
002400 
0C2410 
Oc24u2C 
002430 
O02uu0 
oc24u50 
OC2446C 
OC24u7C 
OC2uAaC 
002490 
0025NC 
0025iC 
002520 
0C2530 
oc2540 
00255C 
002560 
002570 
002580 
002590 
002600 
002610 
Oc2é2C 
00263C 
OC2640 
002650 
002666 
002670 
O0C268C 
002690 
00270C 
002710 





31 
ae. 
33 
34 


340 


35 
50 


51 
52 
53 
54 


SSOCMEGAM(M)=OMEGA(LI+(GMEGA(L-1)-GPEGATL)) ©ABSF(DEL (LI) (A BSF(OEL(L) 
TI+ABSF(DEL(L-1))) 


ce) 
76 


COM=.5*LCM 

L=L+} 

IFC JK) 100C, 1000, 1002 
BM=M 

IF (HM-BM) 3000,3CCC, 22 
M=M+] 

TF (M-3) 38C, 44,34 - 


ADOM=(OMEGAM(N-1)-OMEGAM(P-2)) #. © 


CR=ALCCY/5C, : 
AA=0.0 : 

BR=C 

cc=C 

GG=C.6 

DCM=AC0M 

DELC=0 
OMEGA(L+1T)=CMCGAM(M-1142.*C8 
CMEGAC=OMEGA(L+4+1} 

GO TC 31 

TFC CMEGA(L)-CMEGAG)55,60,6C 


OMEGA(L+T)=CMEGA(LI-COMCCA(L)-CM-CAC 1 /2.C€ 
IFC CMEGA(L)-OMEGA(L4+1)-CRIS1551, 31 


cc=1.¢ 

GO TC 3} 

OELC=VEL(L) 
CMEGAC=OMEGAEL) , 
OMEGA(L+1)=CMEGA(L)+00™ 
GO 10 31 
IFCCELC#OEL(LI129,29455 


GO TO 32 
IF(CMEGA(L)-OMEGAG)21521,22 
IF(BB)58,58,56 


PR=1.0 
TFCOMEGAIL)-OMEGAIL-1)-CR1 25529, 30 
M=M- 1 

GO TO 3600 

IF{GG)70, 70,72 
IFCCEL(LIIT7 71,7? 
GG=1.0 
OMEGA{L+1)=CMEGA(L)-C2 
GC 1C 31 
IF(GG-1.0)73,73,5 7u 
GG=2z.0 


HOEL=CELI(L) 
OMEGA(L+1T)=CMEGA(L) #2. #CR 
co Te 31 

TOEL=OEL(L) *HCEL 

IFC TCEL) 75, 76576 
CREGAM(M)=CMEGA(L-2) 

GO IC 33 

KK=1] 

M=M-] . 


GC TC 20Ccc 


OC272C 
002740 
002740 
0C2750 
B027EC 
002770 
oc272C0 
902790 
O02ECC 
902510 
0c2870 
0C283C 
oo2euc 
002850 
0C2840 
G02E7C 
002e°¢ 
00289C 
062900 
002910 
00292C 
002930 
oc2940 
0¢2950 
00296 
002970 
002980 
0c¢2990 
003600 
003010 
063020 
00303¢ 
003C40 
0C3C50 
003060 
003076 
0c3c80 
003090 
003100 
003110 
03120 
063130 
003140 
003150 
003160 
0c3170 
003180 
003190 
003200 
003210 
003220 
003230 
0C3240 
003250 





AAD 


77 TF(CC)56,56, 54 
OUTPUT 


3000 TFC KMI4CO,4CO,4uCe 
YOO PRINT 4G} 
HOU FCRMAT (IHI//) 
PRINT 4&C2 
HOZOFCRMAT (15H PROGRAM VIPIPCy235Xe37hG Eo FINK NrA-1,359X, 1SHFEBRRLAR 
T¥ 1964//,10x%, 10f&R MOBAL FREQUENCICS OF & PLANAR PIPING 
c¢ SYSTEM ARE LETTRPIRED BY AN ITERATEVE PROCCFEURE USING TRE /ZSH ME 
ZTHCL CF TRANSFER MATPICES.///11GK ae RHR RENEE RRERH HERE ETRE 
WeeTTCTCECECCCCCCCCCCCCOOCOCOT CEOS CeCCCOCEOC OTe Cee ee eee eee ee eee eT 
Cee eee eee eter eee ees 44) 
HO3 TFCNMIUGSE, 4058 ,4uCk 
4O4 PRINT UC5 
4OS FCRMAT (1hF1//) 
4053 NM=j 
WO59 ITFUSKIUCEC,4O40,4CHS 
4C60 PRINT 4UTO Ju 
HOTOOFCRFAT (5CX,8HPROPLEM I2,//75X,47HIN PLANE MODE FRECUENCEES (RADIA 
INS PER SECOND) 4/45 %,4RMOLE, 15%, 10H FREQUENCY /) 
GC TC 4079 
4065 PRINT 4YC7E JS 
4CTSOFORMAT(SCX, SHPROBLEM [2,//32r,S51hCUT OF PLANE MCCE FREQUENCIES (RA 
TDOIANS PER SECCNOC) 4/4USX,4RMCCE, 15%, 10H FRECUENCY 4) 
4O79 MH=M 
PRINT YOR (M,yCMCGAM(M),¥=1,hH) 
HOB FORMAT (46X,12,17X,FIN.E /) 
IF CKK)4C80,4090,408C 
4G&O PRINT 4085 
KOBSOFORMAT(1OX,é69FPRCRLEM TERMINATEC CUE TO SIGNIFICANT FIGURE LIMITAT 
TION Ct CCMPUTER. /) : 
4090 IFCAMYK)412,409,412 
409 PRIAT 412 


4TTIOFORMAT (16X,53H A LUMPED MASS APPROACH WAS EMPLCYED, A 
IND )} 
GO TC 414 

412 PRINT 413 

413 FORMAT (24X,45HA CISTRIBLTEC MASS AP PRCACH WAS EMPLOYED ANDO ) 


44 TF(S0)415,415,420 

KIS FFCRT)416,416,418 

416 PRIKT 417 : : 

BITOFORMAT (2X,67FTFE EFFECTS CF SHEAR JEFLECTICN AND ROTARY INERTIA W 
TERE CONSICEREC. /) 
GO TC &z5 

418 PRINT 419 

HISOFORMAT (3X,63HTHE EFFECTS CF SHEAR DEFLECTICN WERE CONSIDERED WHIL 
YE THOSE OF 43M, Z2S5HROTATIGNAL INERTIA WERE NEGLECTEC. /) 
GO TC 4z5 

420 TFCRT 4C1T 421 sacd 

K27 PRINT 422 

4H220FORFAT (3X,724THE EFFECTS OF ROTATIIJNAL INERTIA WERE CONSIOEREC wH 
TILE THCSE OF SHCAR /3X,CTHDEFLECTION WERE NEGLECTED. ‘) 


D8 


003266 
OC3270 
0C328C 
0035<5C 
Jessa CC 
002210 
Ces sc 
003336 
OC 3540 
003350 
002360 
0033/7C 
003280 
COs 2910 
oc340Cc 
003410 
Oc3420 
0034<0 
OC3K40 
003450 
003460 
003470 
00349C 
003490 
0c35C0 
003510 
00352C 
003530 
003540 
003556 
00356 
003570 
003580 
003590 
003é00 
003610 
0C 3620 
003630 
003640 
003650 
003660 
003€70 
003680 
004690 
003700 
003710 
003720 
0037 30 
003740 
oc375C 
003760 
003770 
003780 
00379C 








=> _ = 





CEC sew 

423 PRINT Uch 

HOUOFCRMAT(E( SK, ,EcETHE FFEECTS OF SRPEA: ULCFLELTION 4NU ROTARY INES TIA we 
TRE NEGLECTEL. /) 

425 PRINT 4uOd 

QHOCKCRMATCIVOH seek ere RRR HERR HR RARER ERE ERR a 
RE ceeCeecrCeeecererceececeee er eereeceereeee en eee cee eee eee se ee eT #eE RRR ) 
TF(GDABC-1.)4511¢452C,44C1) 

QhOT PRINT 4iCe ° 

QHOZOFCRMAT(S5H SECTICN LENGTE ANL SUPSEZ TION WATA/Z2ZXK,14+Stl TICK ALMFE 
TRePAVC SHLENGTE CF SECTICNC INCHES) ¢5X,27HNUMFE® CF SceCTIINS ) 
Beta Ws Cty TLCL ye 24 tbs L=WNS. ) ‘ 

UR3Z FORMAT (29Xe1h 2525XeFT ots clXeFu.O) 


LMEL 
PRINT 4b] 

WHTOFCRMAT (11H GRAPE DATA /2°X, TERT TERATICN NUMBER, 4X, SCHFRECUCNCY (R 
TPADTANS PER SECOND), 3X- 25+ VALUC C& FR. GUENCY CETERMINANT ) 


PRIAT &44 (L,CMCOACLI,CEL(L),L=1,LM) 
WAY FORFAT (29X,13,°4XeF Il lobe 18x,F 15.7) 
"S520 PRINT hy? 
QL2 FORMAT(IX,IGHINPLT CATA s/) 
PRINT 4521 
US2TOFORMATUUXs SFCOCMPCNENT »UxX,SELLAMETER, SX, V4HWALL THICKNESS »5X,6RLEN 
IGTHs5X,9HRALIUS OF o9Xs SHINCLUNEN 5X, (HOENSI VY 5 5X5 7HELAST 10 9 5X_SHSH 
ZEAR /EPX,SHCURVATURE SX, V2FANGLOE CF ARC, VEX5 7RMCOLLUS » 5X y 7 MCEULLS 
B) 
BEC 4525 [=1,NS 
TFOPCTIINS24 bE 2u,uSe 4 
4524 PRINT 4522 TeCl Ls CTCL be TECT) »RHCCID, TEL TROL), BAMUCT) ECT) G(T) 
R5220FORMAT (Ske 13,2 oF Peay IC Xp FE eo 5e Pay Fon Zeb X a FE 5g GX pFO ee Ke Flee euXy 
TEBE. 2,4X, Eber) 
GO TC 4525 
4523 SH=1. 
4525 CCNTINUL 
TF CSh)4S30,4530,4826 
4526 PRINT 4532 
4532 FORMAT (//1X_,7FEANGERS /) 
PRINT 4527 
YSZ27TOFCRMAT (2X, 9HCCMPCNENT 5 TUX, <HELX 9 DYX s ZHCLYV» TUX, ZHCLZ 5 TUX » ZHCTX, TUX 
Te3HCTVe 14 Xe S5HCTZ /) 
4528 CC 4530 [=1,NS 
TFCP CT) V4530,4520,4509 
W529 PRINT 4537 LeCLXCT) eCLY( TV ,CLZ(1) CT XC1),CTY(1L) C1201) 
4531 FORMAT (SXe14_15XeF 10D e TX pF IGe ds MKF VCoc af XeF IOs PXsF IDs 27K oF IC. 
12) 
4530 CONTINUCL 
* PRINT 4he 
4S FORMAT (//2CH BCUNDARY CCACITIGNS ) 
PRINT 4bO (SVICT,1) 1=1,€) 
MUG FCRMAT (5H SVI=6F72.0) 
PQRENT 457 (SVEC(I,1) I=1,€) 
451 FORMAT (SH SVESoF=.0) 
IFONBRIGS11,4511-,45C2 
4502 TF(JK14505,4505,4815 
4505: PRINT &5I00(SVIPC I,J} 1=1,6) J=2sNNNI) 


D-9 


Ceseco 
CC3U1C 
CEB S20 
0c 2% 40 
003540 
N0325C 
003860 
lextomae} 
903880 
0043<¢9C 
0039C0 
0C3910 
aeC3970 
C0395C 
OO39uC 
003950 
0043960 
00397C 
090398C 
002990 
ocado0 
oo“uolc 
OCuCzZO 
0040.0 
OO404u0 
O0k805C 
004960 
Oo4c7o 
oo4oRko 
008090 
oc4l00 
004110 
004120 
004130 
OCkIT40 
004150 
004160 
Ocki70 
0C4 180 
004190 
0o0k200 
004210 
OCKkZ2C 
004230 
004235 
OCc4uz40 
OC4250 
0048260 
004270 
OOwzZaC 
004290 
004400 
OOuZsI1C 
004220 
004340 





4510 


4515 
4517 
u$11 

46 


U47 


452 


453 


FORMAT (¢€F3.0) 
GC To 451) 
PRINT VOITI(SVOP( THs) f=1,0) Jeze KNX) 
FORMAT (46F3.C) 
CONTINUE 

NC Ub6 M=l,hH 
OMECAM()=0.0 
NC YW? L=lylM 
CREGA(LI=C.0 
NEL(L)=0.0 
OELC=C.0 
CMEGAC=0.0 
TFOMND 4S 35,453,452 
MN=O 

L=) 

v=] 

AA20.0 

BB=0.0 

cc=0.0 

GG=0.0 

KK=0 

ADOM=EDOM 

CR=BR 
OMEGA(T)SOMEGAT 
OQMEGAC2=CMEGAT : 
COM=ADOM 

Ge Te 998 
JJzJJ+1 

DO 454 [21,AS 
P(1)=0. 
Z(1)=C.C 

GQ TO 456 

END 


SUBROUTINE SUBSEC (RHO, TFET4 + Ce TL eH4eZ) 
RHOV=ABSF(RHO) 
IF(HM-5.018,8,7 
HM=5.0 

RAZTL/D 
PF(RHOD-0.)3,224 
TF (RA-1.155596 
z=0. 

RETLRN 

IF (RA-3.) 15,175,164 
@=2.0 

RETURN 

IF (RA-6.)17,17,18 
223.0 

RETURK 

TFC HM<3.) 19102 
226.0 

RETURN 

222.0 8HM 


De-10 


004340 


OCW 2506 
004360 
004370 
0C4380 
004390 
oo4uuo0 
Ockutd 
004420 
0044430 
O0uuLO 
o044us0 
004460 
oo4un70 
OC4uyEO 
004490 
o0c4u5co 
004510 
004820 
004530 
ocu540 
004550 
604560 
004570 
OO45f0 
004590 
004é00 
004610 
004620 
0c4ue30 
OOuEnO 
004650 
004660 
ocué7O 
004680 
004690 
004700 
004710 
004720 
004730 
004740 
004750 
OO4uT6O 
ocuT70 
004780 
004790 
004800 
004810 
004820 
004830 
004840 
0048850 
004860 
O04¥ETO 





fFwh 


RETURN 

RN=RHCO/C 
TF(RL-1.)19,19,23 

Z=0 

RETURN 
TFCRL-3.)24,24,25. 
TFCTRETA-.44) 26, 26,27 
z=0. 

RETURN : 
EY SS ei eit ale 5 rely 29) 
2=2.0 

RETURN 

Z=3.0 

RETURN 
IFCRE-6.)30,3C,4 
TFC TFETA-.4u) 22, 
Z=2.C 

RETURN 
TFCTELTA-1.57) 24, 24435 
z=3.0 

RETURN 

Z=4.0 

RETURN 

TFC TRETA-.44) 26, 36,3 
2=2.0 : 
RETURN 

END 


1 
ely oS 


SUBREUTINE DISTM (AMU,ZALY»,OMEGA, IL 2 ,0,019GySDyeRI gE, FRAC pu) 
DIMENSION U(6,6) 
AQA=(3.14159Z2654U/U.0) e#(Cee2-Nleec) 
AREA=ARA®FFAC 
BE=TL#®CMEGA#SCRIF(CAMU/(E#ARE)) 
PU=REAMUeTLe eb ae CMEGAee2 
IF(RT)1,1,2 
T=ReAMUH (ALY # CMC GAeTL ) e¥Z 

GO TC 3 

T=0.0 

TF(SC)4,4,5 

S= (APU (CMEGA#TL) #22)/(G8ARZA) 

GO TC 6 

$=0.0 

SPT=S+T 

VV=SCRTF(PU+(S-T)##2/4u.0) 
AL1=SCRTF(VV-.5*SPT) 


" AL2=SCRTFCVV+.5*SPT) 


CL=1.O/( ALI ¥22+AL24#2) 
CHLI=CCSH(AL1) 
SHLT=SINHCAL1) 
CL2=CCSFALZ) 
SL2=SINFCAL2) 
SPE=SINF(PE) 


CBKE=CCSF(BE) 


Dell 


004820 
OCus8IC 
CO4uSCO 
OCHSIC 
004920 
OGK9ZC 
O0uk940 
Ccugs0 
004940 
OCuSTC 
00493¢C 
004990 
ocsco00 
9o05cic 
0c5020 
005030 
ocsouc 
00505C 
005C60 
005070 
005C80 
005C9C 
0051C0 
005110 
00512C 
005130 
005140 
005150 
005160 
00517C 
005180 
005190 
0c52C0 
005210 
005220 
005230 
005240 
005250 
005260 
005270 
005280 
005290 
005300 
005310 
00532C 
005330 
005340 
005350 
005360 
005370 
005380 
005390 
005400 
005410 





i) 


nO ew 


COO=CLa(ALZeelaeCHL IT +AL Tee? eC?) 


COV=ACLe( CALC ea2eSEL ED /AL T+ Cob tee #5 20/7410) 


CO2=CL#(CRLI-C€CLe ) 
COZ=SCLe(SHLI/ALI-SLE/AL2) 

(i le J=1,6 

OC 7? K=l,¢4 

U(J3,K)=0.F 

LOT, 1I=CrE 

UCT, G)HTL#SPES( CE FARARE) 
U(2,2)=COC-S§*CC7 
U(2,2)=TL*(CC1-SPT*#CO3) 
Ule,b)=ReCOZeTLeae 

Ul], S5)=CReTLee2Z/Pu)e((PY4tSuaz )aCL2-S#CCl ) 
U(2,2)2P4UeCCZ/TL 
U(3,2)=s=CCC-T#CCZ 
U(3,4)=TL#Re(COl-T#CO?) 
U(3,5)=UCZ,4) 
U(R,2)=PURFRROCOL/TL eee 
U(4,5)=Q28Re( (Pht fea2)eCC3-T#COld/ Te 
U(4,4)=U 03, 3) 

Ut, S)=L (2, 3) 

U(5,2)=PhHaRRe (COI-S*CO3) /TL #2 
LI5,35)5L(6,2) 

U(5,4)=U(3,2) ‘ 
U(5,S)=UL02,2) 

Ul6_ 1) =-AMUSTL#CMEGAGH2eSRESBE 
U(6,e)=CRE 

RETURN 

END 


SUBRCUTINE CISTMC CAMU, ATX sATY»,OMEGA,TL,AJT»ReD-D1,G,SD, RIT ,FFAC,L) 


DIMENSICN UE ,6) 

Wwele/(GRAJT) 

ARE A=(3-741592654/4 .O) # (Ce, -Dl**2 )* FFAC 
BE=+SCRIFCAMU # ( TL FAIX#CMEGA) # 2h) 
PUESRHAMU # TL ##haeQMEGAra_ 
TF(RI)2,253 
T=ReAMUS(ATY#COMEGASTL J tee 

GO TC 4& 

T=0.0 

TF($0)5,556 : 
S=(AMUt(CMEGA#TL eat) /(GHARLA) 
co 1c 7 

$20.0 
VV=4SCRTF(P4+(S-T)eec/u.C) 
SPT=S+T 

AL1=+SCRTF(VV-.58SPT) 
AL2=+SCRTF(VV+.58SPT) 
CL=1.O/ (ALT ##2+AL2 #22) 
CHLI=CCSH(ALI) 

SHLV=SEINHEALT) 

CL2=COSF(AL2) 


Dele 


COSueC 
005440 
ocs4so 
OO545C 
OO546C 
005470 
0c5430 
0C54&9C 
0C5£CC 
005510 
005520 
005% 40 
O0osS4c 
6c555C 
005560 
CUSsnc 
005570 
0¢5590 
005600 
ocsé10 
005620 
005640 
005640 
OC565C 
00566C 
005670 
005680 
00569C 
005700 
005710 
00572C 
00573C 
OC574u0 
005750 
005760 
005770 
005780 
005790 
0058C0 
005810 
005820 
005320 
005840 
005850 
OC5860 
0C5870 
0058930 
005890 
0059C0 
005910 
005920 
005930 
005940 
005950 





ew 


SL2=S1Nf (ALC) 
SPE=SINE (PE) 
CPE=CCSF(FE) 
CCO=CL#(AL2¢#28CHL 


oe Wb oo, Dell 24) 


COIACL#( (ALE se2eStLIP/ ALY # (AL Te#2®SLZ)I/AL.) 


CO2=CLe(CRLI-CL?) 


CO3=CL®(SHLI/SALI-SL2/4L2} 


Lo 1 K=1,6 
U(J,K)=C0.C 
L(1,1)=CVE 

UCT s2)=TLeWweSPE/RE 
U(2,1T)=-AsU ® TL 
L(2,2)=CbE 
U(3,3)=CCC-S#C02 


ATX eT eSPE/Si 


U(3,4)= TL #(CC1-Selecets) 


U(3,5)=ReCOzeTLe®? 
U(3,6)=(ReTLee2/Pu 
U(4,3)=PheCOs/TL 
Uth,4)=COC-TFecCe 


pe(-S#CO1l4+(Putesae2 )*#CC2) 


U(K,5)= TLeR* (CO 1-T*C03) 


L(4,6)=0(7,5) 


“ RR=1.C/R 


UIS5s Zd=PHRRRACOZ/TL HZ 


U(S5,4)=RRe(-T#CCi+t 
LI5,5)=U(4,4u) 
U(S,6)=2U(3,4) 


(PH4+Tee2)eCO2ZV/TL 


U(6,2)=PH#RRe(CLlI-S#CN2)/TL 9 #3 


U(6,'4)2U(5, 4) 
L(6,d)=L 04,2) 
U(6,6)=U(3, 2) 
RETURN 

END 


SUBRCLUTINE SFIELC 
DIMENSION A(6,6) 
ARA=(3-34159265u4/b 
AREA=APASFFAC 

OL =TL/éZ 

DO 1 J=2,5 
Al(1,J)=0.0 
A(1,3)=1.6 
Al1l,6)=0L/(C#ARA) 
A(2,1)=6.0 
At(2,z2)=1.0 
A(2,*)=DL 
Al(2,4)=DL ees" R/7 oC 
IF (S0)2.2,3 
A(2,S5)=LLees2#R/é.C 
Ge ic & 
A(2,5)=UL#*#:2R/6.06 
A(2,6)=C.0 
A(3,1)=C.0 


(CT ,TL Py 2Z9GySl 2D,y Ev FAC, A) 


oChe(LeaxZ-Ciees ) 


“UL /(GSARLA) 


0065960 
005970 
00598C 
0059906 
OC6GCO 
006010 
Ococ.C 
0060450 
aceouc 
0C6G30 
aocécéc 
066070 
006080 
00609C 
006100 
00611C 
00612C 
00613C 
006140 
006150 
006160 
004170 
006180 
006190 
006200 
006210 
0062<C 
006220 
006240 
006250 
0062é0 
006.270 
006280 
00629C 
006300 
006310 
0065290 
00633C 
006340 
006350 
0064460 
00637C 
006380 
006290 
006400 
006410 
006420 
006450 
006440 
006450 
006460 
Océ64uTo 
o064eC 
006490 





— 


Ww NO 


A(3,2)=C.G 
AS So) = eo 
AC3,4)=uL#R : 
A(3,5)=A(l,L) 
A(4,6)=0.0 

CG 5 J=l.¢ 
Al(4,J)=C.C 
A(5,J)=G6.0 
A(6,J5)=0.C 
A(u,4)=1.06 
A(4&,5)=hL 
A(5,5)=1.C€ 
A(6,6)=).0 
RETURN 

END 


SUBRCUTINE SFIELC (Cle TL sASI R28 e25G SURI, 0,FFAC,A) 
CIMCNSIGN A(6,6) 
AREA=(3.1U15S926EU/4 CO) aC Lee: -fLer2 je FFAC 
DL=TL/2Z 

DC 2 J=1,6é 

CO «< Kelyé 

A(J.K)=0.0 

A(1,1)=1.0 

A(2.2)21.G ‘ 

4(3,3)=1.0 

A(4&,4)=1.C 

A(5,5)=1.0 

4(6,6)=1.0 

TF(SCORI-O.)7,7,8 
4(3,6)=DL##5#R/6.C-CL/(G#ARLA) 
GQ 1c 7 

AU3Z,E) =f LeeseR/E.C 
A(1,2)=UL/(G#AJT) 

A(3,4)=CL 

AC3,5)=LL¥e2#R/2.C 

A(&,5)=tLeR 

At4,6)=0(3,5) 

A(5,6)=0L 

RETURN 

END 


SUBRCUTING CFIELD (ReZ,6,S0D,RHC,THETA,N,NI,£,A) 
DIMENSTICN A(6,6) 

PHT =THETA/Z 

RO=ABSF{ RHO) 

IF( RHO) T5152 


V=-1.0 
GO Tc 2 
V=1.0 


CP=CCSF( PFI) 
SP=SINF(PFI) 


D-14 


c06£00 


BOOZ LG 


OCb6ac6 
CC6320 
Océtuc 
90633C 
9C6560 
006570 
OCESEO 
006590 
0C066CO 
O0C6610 
0C6620 
O06¢45C 
006640 
9066 30 
006eEese 
O04E70 
006620 
066690 
9067C0 
Oocé710 
006720 
0067 50 
océ740 
006750 
0C6760 


-OC6TTO 


océ7e0 
006790 
océeca 
OO0éEIC 
S06E20 
006830 
O006E40 
006850 
006fE0 
006870 
906880 
006890 
006900 
006910 
006920 
006930 
006940 
006950 
006960 
006970 
0069RC 
006990 
007000 
007010 
oo7Cc2C 
0C7C30 





11 
12 


ee Pe ST ares 

COC 4 K=1,6 

Ndi nS NSE 6 @ 

At(1l,1)=CP 

Atl,2)=-V#SP 

AQ 1,2) =-VeR8ve(1.0-CP) 

AQ), ,4)=-VeRuUseze ce ( PHI-SE) 
AR=S$.VUTE92E Sua (Leat—-CLae7 SKC 
w=T.0/(c#AR) 

FS=(SP-PHI#UP)a.d 
Fl=(PHITeCP4+SP)«#.5 

FSH 22.0-2.08CP-PRET #SP)a.F 
FO=(2OePHIT +PHT eC P—2.0eSP)e.5 © 
ACT SIAVE(IRCthePPlLaSPe. -Ret S82 #5) 
A(T, S)=RLeWwtb 42 eF Oe Qh ae 2 
A(2,17)=V"SP 

A(2,2)=CP 

A(2,3)=hD«SP 
A(2,4)=ReERDaeZe(1.0-CP) 
A(2,5)=R20eF “eC ntReRCaeZ) 
Al2s6)=All, 5) 
A(32,/)=1.0 
AU3S,4)=ROCaR # PEI 
A(3,5)=A(5,4%) 
AL3Z,E)=A(1,4) 
A(4,4)=1.9 
A(4U,5)=H4(2,7) 
A(4,6)=A(1,2) 
A(5,5)=UP 
Aloe niles c 
Al6,5)=A02,1) 
Al(6s6)=CP 
RETURN 

END 


SUBRCUTIAE CFLELC (AST:9nZ2eCySCRIyRHC, THETA, DCI, FF AC, A) 
DIMENSICN Al6é,5) 

PRTI=THETA/S2Z 

RHCLSREC 

TFCRHC-C.O)1ts1),12 

V=-1.C 

Ge TC 2 

V=1.C€ 

CP=CCSFCPRT) 

SP=SINF( PHI) 

w=1.0/( G #AJT) 
FI=(PHIeCP+SP)/2.C 
F3=(SP-CPePHI)/ 2.0 
FS=(2.0-2.08#CP-PRI*SP)/2.C 
F6=(2.0aPHI+PET&CP-?.C8SP)/2.6 
RHO=ARSF (RHC) 

CO < J=1,6 

DO z K=196 


Del15 


Og7C 4 
OMe Ble, 
co7ceg 
O0TL70 
a) 77 @ (C 
oc7c9c 
oo71ico 
9c7.10 
otal ete 
OST Zl 
CO7140 
OG files 
OC7iec 
o07TI7C 
co7ica 
ole Wa mele 
Oc7.2.C0 
Vcr es 10 
OC maa 
007.450 
OCTZ40 
GT A S310) 
007260 
OOO 
oc7230 
OC7290 
Q0074C0 


OOTSIG 


OC VeseaG 
007230 
OC? <40 
oo7ssC¢ 
OCc7260 
007270 
007480 
oc7s9c 
ocT4o0 
007410 
OC7THZC 
907450 
oOcTa40 
007450 
OCTHEC 
007470 
oc74uad 
007490 
007f00 
007510 
ac75°0 
Ooc7T53C 
O07TSu0 
007550 
067560 
007470 





14 
5 


-Wwh 


A(JsK)=C.6 

IF (s0R0)9127,15,14 
APEA=(3.1415926S54u/4 .C)e(Ce# -BCian2 je FE AC 
SN=RHC#THETA/ (GaZe#AREA) 

ti te aS 

Sb=0.C 

A(1,1)=CP 
A(1,2)=(wWeRHGsF1)-(RHOeF ZeR) 

A( 1,4) =VeSP 

Alls, S)=Vel(wteR) #(REOePHI®SP1/2.0 
A(T, 6) =Ve(WeR) o( areas) ar? 
A(2,2)=UP 

A(2,5)=vesP 
A(2,6)=V#RHC#(1.CG-CP) 
A(3,1)=-A(2,6) 
A(3,2)=-1.08A(1,6) 

4(3,5)=1.C 

A(3,u)=RHC#SP 
A(3S,5)=(CR#(REO##2) aPHE®SP)/2.C)-WeR Cee cer 
A(3,6)=(R#RHO4e 72 #F 3)-(WeFeeRbCee 7)-SN 
A(4,1)=-V#SP 
A(&U,2)=-Vel(weRPeREOePHI @SP/ 2.0 
A(&u,&)=CP : 

Ath, S)=(RARHG AFT) -(WeREC FS) 
A(&,6)=4(2,5) 

A(5,2)=-V#SP 

A(5,5)=CP 

A(5,6)=A(3,4) 

A(6,6)=1.C6 

REO=RHCC 

RETURN 

END 


SUBROUTINE PCINI (AMU,ATY,QMEGA,TL»Z,R1,8) 
OIMENSTON 6(6,6) 
AM=AMURTL/( 2-1.) 

LC 1 J=1,6 

DC 1 K=1,6 

B(J,K)=0.0 

B(1,1)=1.0 

P(2,2)=1.0 

B(3-5)=1.€ 

TF(R1)2,2,3 

Bib ,2)=-AM# (ALY @CMEGA) #82 
B(4,4)=1.0 

6(5,2) =AM*OMEGA# eZ 
B(5,5)=1.0 

P(6,1)=-B(5.2) 

8(6,¢6)=1.0 

RETURN 

tND 


D=16 


‘OC7526 


oc7réegc 
oo7éece 
CO7e1C 
OC7h=C 
CO7¢ SC 
OoTEe4NC 
C7650 
Ooc7é460 
oo7¢e¢7C 
OOTéELC 
C0769C 
oo7v7cc 
oc77ic 
oo772C 
OC isc 
QOTTRC 
Cc7750 
007760 
oc7770 
007780 
ooT7S0 
o078eo0cd 
0c7810 
007820 
oc7T#350 
oa7exud 
oc78s50 
007860 
OO7E?O 
007890 
oorTeg9a 
007900 
00791¢C 
007920 
007930 
oc794u0 
007930 
0c7960 
oc7970 
oc7920 
007990 
ooscco 
0c8scl10 
oo08c20 
0080350 
008040 
008050 
008060 
008070 
ocsosc 
008090 
0c8100 
008110 





SUBROUTINE PCINC (AMU,AIX,ALYs OM: GAy Thy ZySDR1yP) 
DIMENSION £&(6,6} 
AM=AMURTL/(Z-1.0) 

DC 1) J=i,6 

DOD 1 K=lyé 

Rit) =Cc 

R(1,1)=17.C 
B(2,1)=-AMB(ALKR®OMECA) eee 
B(2,2)=1.0 e 
F(3,2)=1.6 

B(4,4)=1.0 ; 
LF(SORI-0.)9,9,1G 

R15 ,4)=-AMR(ALVACNEGA) #82 
GQ TC 11 

b(5,4)=0 

le (So Sey oe 
P(6,3)=AMROME CAN eC 
P(6,6)=1.0 

RETURN 

END 


SUBROUTINE MATMUL CAsB,Z,U) 
DIMENSTON A(616),F(6,6)20(616) 501695) 2U( 656) 
be 1) J=1,6é : : 

LO 1 K=1,é 

C(JS,K)=0.C 

bO 1 L=1:6 

C(S,KIV=C(5,K) 480 U,L)#A(L,K) 
N=Z-1.8 

(Dla) 2 J=1,6 

DO Z K=1,6 

U(J,K)=C(J,K) 
IF(Z-2.0)6,6, 2 

DC 5 M=1,N 

CO &Y J=196 

OO 4 K=1,6 

D(J,K)=0.0 

CO 4 L=1e6 

DUI pKIHLES,KI+C(SeLI#UCL eK) 
O00 5 J=1,6 

CC 5 K=)e6 

UCJ,K)=005,K) 

DO 7 J=126 

OC 7 K=1,é 

O(J,K)=0.€ 


' 00 7 L=l,6 


COIeKV=HL CI K) FAC Sy L FUL, K) 
OC y J=1T,é 

DO & K=176 

UlJsK)=O05,K) 

RETURN 

END 


Del? 


ocel20 
C0813C 
OCckixe 
0C3150 
OcBi¢d 
OV 
0c81iz20 
0cel90 
OC8=ECO 
0C8z2 10 
oc8220 
0C8&<z 40 
008.40 
0C8250 
008260 
008270 
OB. 2C 
008z9c 
oce<s00 
008216 
008420 
008330 
008240 
008350 
008460 
0084570 
OC8380 
008390 
008400 
008410 
068420 
008430 
008440 
0c8450 
008460 
OC84u70 
0084 2C 
008490 
008500 
008510 
008520 
008530 
008540 
00855¢ 
008560 
008570 
008580 
008590 
008600 
008610 
008620 
008é3C 
008640 
008650 





C 


Fw 


30 
31 


SUBROUTINE RICIE CANL, TL»CMEGA,AIYy2%1,L) 
DIMENSION U(656) 

AN=AMURTL 

nc t J=l,é 

CC 1 K=1,6 

(i (Cal g (5 SUE RICE 

U(1,1)=1.0 

UE, TL) =-AMeCYEGARR?- 
U(6,6)=1.C 

U(2.c)=1.0 

Ul2,5)=TL 

U(3,3)=1.C6 

UCU 2) =AMaTLeCMEGAaec s/s? JC 
TF(R1)2,2,3 

UE, Ss) =AMaOMECAR ace (TL ent /€.C-Alye#? ) 
GC TC 4 
U(4,7)=AMe(CMEGARTL I 8 82/€.C 
L(4,4 )=1.0 

U(4,S5)=IL 

Uu(5,2)=-U(6,1) 
U05,3)5b(4,2) 

U(545)=1.0 

RETURN 

LND : ' 


SUBRCUTINE RIGIC (AMU,TL,AIX,CMEGA,ATY,AJT,SUR1,U) 
CIMLENSION U(6,6) 

oe . JE ao 

£0 2 K=1,6 

U(J»sK)=C.06 

AM=AMUETL 

U(1,1)=1.0 

U(2, 1.) =-AMUETLa( AIX #UMESA) #42 
U(2.2)=1.0 

Uu(3,2)=1.0 

U(3s4)=TL 

U(4,u)=1.0 

U(S,Z=AMeTLeCMEGAS#2/72.C 

IF(SURI-C.) 36,70,31 
U(5,G)=AMeQOMECAs aZa(T Lee c/6.0-ALY*##2 ) 
cc ic 1 : 
U(S,b)=AMeOMEGAsele(TL#*#2/46.0) 
U(5,5)}=1.C 

UI5,6)=TL 


U6, SI=AMHCHEGAaRE 


U(6s4)=U(5, 3) 
U(6s6)=1.0 
RETURN 

END 


SUBRCUTINE STIFCC (THETA,RECsU) 


D-18 


908660 
0C8570 
008620 
O08e9Cc 
9CErCC 
acs7ic 
008720 
OUg\ss IE 
OC874C 
0¢875¢ 
OCETeC 
008770 
0087sa 
ocs790 
Ocgsca 
O0BY1C 
908820 
nC8e270 
Oc3c40 
aceese 
ocesse 
ocse7o0 
ocsEeEec 
008e9Cc 
ocs9cc 
008910 
acg9-c 
008940 
008340 
008950 
aceséc 
oces7c 
OC896C 
008990 
00900C 
009010 
0c9go2Cc 
009C43C 
OOS0LC 
0Cc9050 
0090450 
0c9070 
o0o9cec 
acgycgc 
9091700 
oc9giic 
00912C 
009130 
0091740 
0c915C6 
00916C 
0c9170 
o0c9180 
0c9190 





AoA 


Po) 


Ww iN 


DIMENSION U16,6) 


(erg J=tyt 
[Ese en K=1,6 
U(JdsKI=0.C 
ACS Ne G4 eS 
Ve-1.0 

Ge loa 
v=1.0 


CT=CCSF(TRCTA) 
ST=SINF(TEFETA) 
U(r, ))=CT 
Utlh,2)=-VesST 
Ul2,1))=vesT 
L(2,Z)=CT 
Wigs = 16 
U(4#,4u)=1.C0 
LU(5,5)=LT 
U(5,6)=-veST 
U(E,f)=VeST 
U(6,E)=CT 
RETURN 

ENC 


SUBROUTINE STIFCO (THETA,REC,U) 
CIMCASICN Ul6,6) 
[TF CRRC)1,1,2 
=-1.95 
GO TC 4 
V=1.0 
CT=CCSF(THETA) 
ST=SINF( THETA) 
OO & J=l1eé 
CO 4 K=1,6 
U(J»K)=0.6 
U(1,1)=CT 
U(1,&)=STaV 
U(2,2)=CT 
Ul2,5).=STe&Vv 
U(3,%)=1.C 
U(4, 1) =-ST#V 
U(H,4d=CT 
U(5,2)=-ST#v 
U(S,5)=CT 
U(6,6)=1.0 
RETURN 
END 


SUBRCUTINE HANGER (CLXeCLYsUTZsL) 


CIMENSION Ul6,6) 
OO 1 J=i,¢é 


CG 1) K=1,¢ 


D=-19 


9cg9.00 
aca. ic 
OC92.C 
4a Ske) 
909240 
OO). Se 
ORSSEC 
OC9L/C 
OCSE0 
0092790 
MOHD ACI 
(overs) = Bf) 
009 °20 
009.76 
0C9354C 
OCIS 0 
o0c9-60 
CCI eS 
009580 
0c9:90 
oc9400 
O0941C 
OCc94zC 
0094 40 
ocouyLC 
OC9YuSC 
009400 


009470 


0c94uc>o 
909490 
0095C0 
009510 
0c9+20 
909540 
OG9I540 
009550 
ocgtéc 
00957C 
009530 
009590 
009éCC 
009010 
009670 
009440 
acgé4uc 
009650 
009660 
009670 
OO9CED 
0C969C 
0c9700 
009710 
o0g7zCc 
O00973C 








— 


U(d,K)=0.0 9C9740 


NA eae 009750 
Cie c= 00 009760 
Wisi V= 1.0 0c9779 
U(4, 3}=CT2 ocgrec 
UE, 4}=1.0 0c979¢ 
U(5,2)=-CLY oc9sco 
US.5)=1.C 009510 
UIG, T)=CLX aegcc 
W16,6)=1.0 009830 
RETURN 0C9£40 
END OCIEEC 
O09YEC 

00927C 

SUBICLTINE HANGEC (CTxXsChLZeCTY el) 009230 
DIMENSICN U(616) 00989C 
ie Veda 15 é 00990C 
CC 1 K=1e6 209S1C 
LIS, K)=0.6 009920 
aie Wil C 009950 
U2, 1I=CTX 00994C 
WZ ecy=1.0 0c9950 
tio 1=1.0 969950 
Ube d=1.6 009970 
U(S,4)=LTY 0099°0 
U(5,5)=1.0 009990 
U16,5)5-CLZ 010cc0 
L(6,6)=1.0 010cICc 
RLTURN 01¢¢20 
END 010C30 
010C40 

010650 

SUPRCUTENE PRANCH (SVE sVVyPEIyMIsU) 010060 
ODIMEASIGNSVE (Se 22) e VVC Gr Eb aU (Or 6 De RVC 25505) hl (25) 9 MM (25 De R202 93) 010070 
Sys ieGlt2.315C20452) 010080 
vel 010¢SC 
DC 2 N=196 010100 
IF(SVB(N@MI)) 25459) 010110 
Ome =1,¢ 01012C 
VVCLeMI=VVEL OND 010130 
=M+] 010180 
CCNTINUc 010150 
DO 4 L=1,3 : 010160 
RU V=VE ly 1) 010170 
RIC2sLI=VVCZyL) 01018C 
RV¢3,Ld=vVC2,L) 010190 
R2¢1eLI=VVC4yL) 010200 
R202, LI=VV(5.L) 010210 
R2(S,L)=VV(E_L) 010220 
CALL INVERT (R1,2eD,LLy YM) 010230 
CCS J=123 010240 
Co > «=1,3 010256 
R(JeK)=C.0 010260 


Oe 3 iheilos 016270 


D-20 





5 


10 


ROJsKIV=El (Se KARO CS LI AQVIL,K) 
SP=SINF( PHT] ) 

CP=CCSF( PHI) 

Cl(t,1)=cP 

G1(1,2)=-Sp 

Gi Gils SpSse. 

Sle,1)=SP 

Glleyi)=CP 


Gi(é,4)=0. 
GU Se pee. 
C1(3,2)=0, 
Giles |= le 
G2(1,1)=1, 
G2(1,2)=6. 
G2(1,3)=0. 
G2(2,1)=C,. 
G2(2,2)=CP 
G2(2,%)=-SP 
G2(+,1)}=0. 
G2(2,2)=Sp 
G2(5,3)=CP 

00 ¢ J=1,? 

DO ¢ K=1,3 
R1(3sK)=0.0 
PQ 6 L=1,3 , 
RUCSeKIERI( Se KI +C2( del dE R(L,K) 
00 7 J=i,3 

He ¢ K=1,2 
R(J»K)=C.0 

DO 7 L=1,2 
R(UJsKIPER CIs KIFRI(J,L)#CT(L,K) 
OO £€ J=1,¢ 

00 & K=l,é 
U(J,K)=C.C€ 
U(T,1)=1.6 
U(2,2)= 0 
U(3,4)=1.C 
U(%,4)=1.0 
U(5,5)=1.C 
U(6,6)=1.0 

DO 9 L=1,3 


U(4,LI=ER(1,L) 
U(5,L)=R(2,L) 
U(6,LI=R(32,L) 
RETURN 

END 


SUBROUTINE STAVEC(SV,N, BC) 
DIMCASICN BC(6,22) 

CG 10 J=1,6 

BC(JUyN)=C. 

IF(SV-1.)151,2 


BC(4,N)=1.0 


De-21 


STC 220 


010290 
01710400 
D101 
@Ge 2 
0103240 
O1C44C 
01035C 
C1076¢6 
010370 
O104#0 
010390 
010460 
010410 
010470 
010420 
C10440 
010450 
O104éC 
0104870 
O10480O 
010490 
01050¢C 
910510 
010520 
010530 
010546 
010550 
9105606 
010570 
0105°¢ 
010590 
010660 
010610 
010620 
01710430 
010640 
010650 
010660 
010670 
010680 
010690 
01707c0O 
010710 
010720 
010730 
O1074C 
010750 
0107¢0 
010770 
01078C 
010790 
010800 
01081¢c 








ww RD 


C= ale) 
EC(E,NI=T.0 

AL TURK 

IF CSV-2.) 4, Sah 
FC(C1,%s)=1.0 
ECtle,h\)=1.0 
PCa yNd=1.0 
RLTLRA 
LF(SV-2.)57596 
PC(s,N)=1.C 
PClusN)=).0 
EC(o,K)=1.0 
RETURN 

IF (Sv-4.d 7, 7k 
RCC1TsNI=T.0 
BC(2,N)=17.0 
RC(L»N)=1.0 
RETURN 
PC{2,N)=1.0 

BC (:.N)=1.0 
RC(6,K)=1.0 
RCTURK 

CNO 


SUBRCUTISE ERANCC (SVP svVV,PEFI,¥1,U) 
ODIMENSIONSVE (Oe 22) US Cd eVV(G Ode WE Ge 41RD eb ec SD yR20 3, FD rRI 3 e5 ds 
GU (Ser) G20 5,7) SC3eZV LEC 25) eMM(25),SRC6,220)eLVI Ee 9) 


v=] 

SROieMID=SVA (191) 
SR(2eM1)=SVe(2,"T) 
SRCZeMTV=SveCh YT) 
SR(4 MI =SVECZe1) 
SRS sMIV=SVEC5,"1) 
SRESeML)=Svel(EeM I) 
UV(O1,1)=vvVC1e1) 
UV(1,2)=VV01, 3) 
UV(1,22=VV(1.4) 
UVE1T»4I=VVE 1,2) 
UV(1,3)=VV(1,5) 
UV(1,6)=VV(1,6) 
UV(.,1)=VV( 3,1) 
UVle,l)=VVI 2,2) 
UVOl,4)=VV( 3,4) 
UV(2,4)=VVI5_2) 
UV(o,5)=VV032,5) 
UV(2,6)=VV(2,6) 
uUv(2,1)=vv(4,1) 
UV(3.2)=VV(4, 2) 
UV(2,7)=VV(4,4) 
UV(3.4)=VV(4,2) 
UV(5.5)=Vvv(4,5) 
UVESeEVEVV(4,E) 
UV(4,)I=VV0e, 1) 


De22 


fal] OVS) 218) 
CHCAsc 
O1dsuc 
010839 
0104450 
O10L7¢ 
(0) | (Ne 200) 
910290 
010900 
910970 
01052¢ 
Cieee Ge 
orcguo 
OVOosC 
910960 
010979 
010920 
01099C 
011CCO 
o11c19 
Oeze 
0110450 
OTL 
011050 
c11ceéo 
cC11c70 
011030 
011090 
011100 
011110 
011120 
011130 
O11)4C 
011150 
011160 
011170 
O111&C 
01119¢ 
011200 
01121C 
Olle ¢ 
011230 
017240 
011250 
O1l2éC 
011270 
011250 
01129C 
011300 
01121C 
011320 
0113520 
011340 
0117350 





ho 


UVt4,2)=VVE2,2) 
UV(4, 3)=VV(2,6) 
UVa ,4)=VV(l,2) 
UV(4_: bFVV0e55) 
UvV(4,é)=VVO_,¢) 
UV(>e1I=VVb5_ 1) 
UV(s,2)=VVC5, 2 
UV(o,.J=VV0554) 
Uv(S»eI=vV(o,o) 
UVOS S)=VV(2,5) 
UV(i,6)=VV05,6) 
UV(e,lIevve, 1) 
Uv(trfd=vv(e, 2) 
UVle,7)=VV(6,4) 
UV(Es4)=AVV(6_2) 
UV(4,5) =Vvlé,5) 
UVEEGEIEVVO EE, 6) 
NC os N=lyé 


LTFCSR(N MEI) 4—351 


(iia! A L=1,¢é 
UVILeMI=FUVEIL,A} 
M=M¢1 

CONTINUE 

00 4 K=T,2 
RIC1sKI=UVO1,K) 
RieeK)=UV6 Ed K) 
R1(5,K)=UVO5,K) 
R201,K)=UVO4,K) 
R2(ZyK)=UV(5_K) 
R2(s5sK)=UV(E,K) 


CALL INVERT (R21,3,0eLL,hM) 


OC s J=1e3 
CC 2 K=1,3 
R(JekK1=C.5 
DC 5 L=195 


RO JeKIER( Se KPERZ(S,LI#RI(LK) 


SP=SIENF(PHE) 
CP=CCSF(PHI) 
Gi(1,1)=CP 
G101,2)=0. 
G101,3)=-SP 
Glle,1)=C. 
Gllcyz)=1. 
GI(2,3)=0. 
G1(3,1)=Se 
G1(32,2)=0. 
G1(3,3)=CP 
G2(1,1)=CP 
C2(1,2)=SP 
G2(1,3)=0. 
G202_,1)=-SP 
G2laez2)=CP 
G2(2,3)=0. 
G2(4,1)=9. 


V114456 
NES 
mal ita 
V1TS3¢ 
C114C9 
CHARS 
O11 ee.S 
O11 peso 
0114uG 
Oe sS 
9174650 
O11W708 
911e7Cc 
Q11e90 
411500 
O11310 
011520 
911530 
011540 
O1155 

0141260 
Clls7o 
011580 
9111S 90 
9116CO 
011610 
OlMeac 
017630 
013649 
011650 
C11éE0 
O11l67C 
O1108c 
011690 
011700 
o1171iC 
ol/rcc 
O111?20 
011740 
COMMAS 
O1l?oO 
011770 
ort7e80 
01179¢C 
oleco 
011¢10 
011820 
011830 
011840 
011850 
011260 
011870 
0118°0 
011290 





10 


wi ho 


G2 :92)=4. 
GZS52 )=)- 
fe « = leo 
Ce a velas 
Ri (upkK)=C.0 
COE = 14.5 


QPOs KIERT OS KItCE(S,LIHRIL,K) 


Wie ¢ Seis Ss 
flr’ if PEP Gg A 
S(JekK=U.C 


S(SeKI=ASCS KVR UC IS LIeTI IL Xd 


sien ao) 


SUPRCUTINE STAVEC 
DIMENSION BC(6,22) 


DO 10 L=13;6 
EC(L»NT=C.0 


IF(SV-1.0) 1,192 


EC(2,N)=1.0 
BCls,NI=1.0 
BC(6,N)=1.0 
RETURN 


TF(SV-2.0)3,344 


BC(1,N)=1.0 
PE (o,N)=17.0 
bC(4,NjJ=1.0 
RETURN 


TF (sV-3.015,596 


BC(Ly»N)=1.0 
BC(4&,N)=1.0 
BC(6,N)=1.0 
RETURN 


IF CSV-4.C)757)5 


011900 
911916 
G119°C¢ 
011930 
011940 
911959 
COMP ee 
011970 
011920 
c11990 
OIZEES 
012C10 
C12620 
012C30 
012CuC 
012c5¢ 
012c6é0 
012C70 
OAC 
012C9¢ 
012100 
012110 
012120 
012149 
01214¢ 
0712150 
012160 
012170 
012180 
012190 
Czas 
012210 
012220 
012230 
012240 
012250 
012260 
0122 7C 
012280 
012290 
0122C0 
012451C 
0312320 
012550 
012340 
012350 
012360 
012370 
o12238C 
012390 
0124uCC 
012410 
912420 
012436 








em 


i 


10 


45 


Gt ta Se 
PCS BSUS 1 
(Clo, d=} 
PAE a Aes 
ENC 


‘ 
at 
20 
" 
av 


SUBIEV TING Ti VEST Abii bok) 


PROGRAM FOR FINCIAG TEL INV RSC °F A RXW 


DIMCASION MB s5 25 Vail Ce) alt. 3) 
SlEAGGis JPRS (LANSKEE SU ILE WIE 
p=1.U 

MORO KETAN 

LOK) =" 

M(K) =k 

BIGu=A(kKyk) 

(CES [ak yh 

LP2G Jekyh 

TE CAFSFOCC IGA) -APSF(ACT,JI1)2 10,2.,29 
PIGA=A(1T,J) 

L«(K)=1 

M(K) =u 

CANT INUL 

INT + RCHANGE 2&CWS' 

J=LCK) 

HF ULURSE 1) S85S80e8 

CO3G [=]l9f. 

HCLL=-aA(K,1) 

ACK, T)=8J,1) 

AC J,T)=t CLD 

INTERCHANGE CCLUNAS 

T=M(K) 

IF(M(K)-K) 45,45, 77 

CNG J=leh 

HOLE =-AC I,K) 

ACSyK)=A 0 Se 1) 

A(J,T)=HOLD 

DIVIttE COLUMN BY FINUS PEVCT 
COS> [=1,N 

Pe i— yee, ooy SO 

ACL, KI=H4ACT, KI /(-A€K,K)) 
CONTINGG ; 
REDUCE MATRIX 

LO6> [=1,! 

C6. vEls™ 

UPL Hats) Se yi Gta 

IF (J-K) 60,€5,60 

ACT, JI=HACT KP RACK, JI+4¢00 2) 
CENTIARUL 

OIVICE ROW BY PEVCT 

CC7h Jetyn 

LECJ=Ky GG, 05 7C 


ACK, J)=A0(K,U) SAL K,K) 


De25 


MATR2I4 


Ozuue 
Ol245c 
21244C 
912479 
012870 
GV249C 
Was Se) 
12570 
Gy Lec 
012°3C 
91254C 
DY25R@ 
Onl s.0 
al 25) 116 
OW 25 
912590 
0126C9 
CZeomllO 
C2226 
0124629 
012¢40 
ON 2e 50) 
012660 
012070 
O22 
012070 
Onane 
OVA7 TC 
Oni aaC 
312740 
G12740 
O27 20) 
012760 
912770 
012780 
012790 
O1l2¢cc 
012310 
01282 

012230 
O12: 40a 
O}es= 

012640 
Cizaan0 
(REE 
012990 
0129C0 
01291C 
9129626 
012920 
012940 
01295C 
012960 
012970 





oy 
ws 


62 
790 


al 


71 


Ch 0 hin 

CONTINUC! PRC ECUCT OF PIVOTS 
D=b4A(K,K) 

RT PLACE P1¥CT RY 
SHOALS Res WLC 9) 
CONTINU: 


Pet ee iL 


FINAL 2%ch ANE COLUMS INTEPUFANCE 


K=N 

K=(*-1} 

WG Ui) YS Ges USitirs tee 
[=L(h} 

TF C1-K) 120,120,105 
OCC J=1,N 

HOLL HAs Kk) 
ACJyK)=-alJ,1) 

ACJ, TY=HCOLD 

J=M(K) 

fFCdJ-K) 160,100,125 
LOlSG L=1,N 
HOLL=AACK , 1) 
A(K,1)=-A( 5,1) 
ACJ,1T)=FOLO 

60 TC Ki 

PETURN 

LAL i 


SUFRCUTISE FEINMAT (Usgtb,I,Vv) 


DIMENS TON Ulle,6),0U lc, 6) ,V(6,6) 


IF CT-V)21,6 26) 
LC 7C J=1,4 
PEO K= 1553 
ViJsKI=U0U,K) 
RETURN 

CO 71 J=l,e¢ 

Ce fl K=1,6 
ViJskKI=C.C 

CO 71 L=1,6 
VOIsKI=aN(S,K) 4UC J, LI eUUIL,K) 
RCTURN 

END 


SUBROUTINE FINBRA (Us VV, RARER V) 
DIMENSION UlG se} ,VV06r 6) VI 694) 


[FENNPR-1)15153 
DORE J=1,6¢ 

LG 2 K=1,€é 
VOJ,K FUCK) 
RETURN 

LO 4 J=1,¢ 

CQ 4 K=1,6 
VOdsK}=0.C 

C9 4 L=l,¢ 


omic 
12970 
SWIC 59 
CARING 
CHUSiGEC 
SU Se) 
913040 
O1305C 
C13¢40 
| Shu, 0) 
0145046 
21507C 
G131Cc0 
ouside 
One 
Culesies.© 
ONG 
Cisse 
YS ee 
OS TieaG 
013134¢ 
012199 
013200 
014. 1C 
Mise 2 
014220 
013240 
Os, 39) 
013.60 
Osis 70) 
Oi sase 
013290 
013500 
013310 
ONE ac 
015430 
012340 
013350 
(a) | Shaker} 
013379 
0123290 
012390 
013400 
013410 
013470 
013440 
G13u4C 
01345¢ 
G13460 
014476 
C1Z24°%C 
013490 
012500 
014510 





Sal 
90 


92 


VU SsKIEVO Se KIFU (UL D#VV(L,K) 
RETURA 
END 


SUBRCUTINE CELMA(LUsSVI,SVEsV) 
OCIMENSTICN UUlEs 6) VV 065615 V (646) SVE CEs 1) SVE(Es 1) 
M=Q 

De 90 N=1,6 : 
IFCSVI(N;1))1,96,1 

M=M+] 

NC 91 J=1,6 

VVC JeMI=EUU( LN) 

CONT INUL 

L=6 

ie 9? K=1,6 

HEU SME CR 5 UNG 52 aa 

L=L4¢1 

De Ge tsi 5 ot 

VIL sKI=VVIN,K) 

CONTINUE 

RETURA 

END 

END 
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013520 
OMS Sesi0 
013540 
C14550 
013560 
013£70 
ONBSeC 
0145590 
013600 
013¢1C 
013620 
013630 
013640 
013650 
0136éC 
013670 
013620 
013€90 
013700 
013710 
013720 
015720 
C13740 
013750 





APPENDIA & 


ACCURACY OF METHOD 


Eel General 

Appendix E contains the vibration analyses of the 
several typical piping configurations to which reference 
is made in Chapter 3. The analyses are grouped according 
to the ourpose for which they were performed. Yhe firs; 
group was performed to establish the accuracy of the method 
embodied in the program VIPIPE, the second and third to 
establish its integrity. 


E.2 Method Accuracy Analyses 
i.21 Systems Analyzed (straight sections) 


System 


Parameters 4 7 ae 3. 
end conditions Tixed-fixed 


fixed-free | fixed-vropped 
length (in.) 
diameter (in.) 


all thickness (in.) 


density (1b./cu.ft.) 


shear modulus (psi) 


elastic modulus (psi) 











weee Cosnarison of Natural ivenucrenees 


@o System 1 


frequenc rad./sec. ) 


Mode 
VIPIPE Comparison Values Difference (%) 
number 




















1 4853300867 48 533126 .00024 

a 133.78346278 133.78345 200000 

3 262 .26843097 262 . 26886 2000016 

4 433. 54478122 433 ..54393 

5 647 .62290622 647 .63991 

bo System 2 
frequency (rad./sec.) 
Comparison Values 

nugaber) 


















7 62356485 
47 679657598 
1338363688. 
26226528107 


iL 
2 
3 
by 
5 433 © 54435149 


Co system 3 


Mode frequency (rad./sec. 
VIPIPE Comparison Values Difference (%) 
number 


7 .6270709 
4706797761 
13383635 
262 26544 
433 0 54415 


204595 
000247 





















°90001 
200006 





000004 








1 33 A4525745 33 445600 000102 
(2 108 . 38602303 108 .38599 000003 
3 22098387715 220 etaea0 200000 
4 386 ..71069206 386.71057 00003 
5 Se os eras eysial 590.10164 ©00012 


* 5). of comparison frequencies for probleas 1 
through 3 are given in E.23. 
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E.e23 Source of Convarison Frequencies for Systeus 1 through 





The governing fourth order equation of a homoreneous 
Single conponent straight section systen 


2 
yiv a Ae y 
Bd 


was solved to give 
y = a-cosXxx + b-coshaxx + c-sinkx + d+sinhaxx 


Upon applying boundary conditions for problem 1 
(y=y! =0 at x=0, and x=L) and simplifying it was found 
that the eigenvalues of the system must satisfy 


cose coshx =1 (Ea237 LD 


2 
soe | 22 Te a 
BJ 


Upon applying the boundary conditions of problem 2 


where 


(y=y'=0 at x=0,y"=y"=o at x=L) it was found that the 


eigenvalues of the system must satisfy 
cosxXcoshe =-1 (B.23=2) 


Upon applying the boundary conditions of vroblem 3 
(y=y'=0 at x=0 ,fy'tO at x=L) it was found that the 
eigenvalues of the systen must satisfy. 
sineccoshet <-cosee sinhoc =0 (2.23—3) 
The roots of equations E.23-1, E.23-2 and E.23-3, and 
from them the first five natural frequencies of the systems, 


were obtained through the use of a digital computer. 





Ho24 System Analyzed (section of constant curvature) 


System parameterss  _ 
radius of curvature 






included angle of arc 
diameter (in.) 

wall thickness (in.) 

elastic modulus 

shear modulus 

density (lbs./cu.-ft.) 

end conditions: fixed-fixed 


Results: - 
in-plane vibrations 












Mode frequency (rad./sec.) 


VIET Comparison Values 
number 
73232973072 73.06 


out-of-plane vibrations 


(rad./sec.) 
Comparison Values Difference (%) 






Difference (/) 








|e 


eo Method Integrity Analyses 


The first group of analyses considered here constitutes 


a mirror image check on system integrity. Four configura- 
tions were considered, and several mathematical models 


were used. The second group consists of two analyses of 


This value was obtained using a curve in reference 2 
which was constructed through the use of a modified Rayleigh 
method. In the words of the reference; "...laboratory tests 
have shown the curves to be substantially correct." 

Comparison values obtained by employing the equation 

W =F ES, after obtaining F from curves in ref. 2. 
uk 
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the same single branch system wherein first one member then 
another is considered to be the branch. 


Eo31 Mirror Image Analyses 
&@ System 1 


system 


Sopn5 
THY Ss NN 


mirror image 


System parameters: See system 2 E.21 


Mathematical model: Distributed mass with the 


effects of shear deflection and rotational inertia considered. 


Results: 


frequency (rad./sec.) 
Mirror Image Difference (%) 


7 .622895556 762289556 200000 


47 .76681871 47.76681872 ~00000 
13363614681 133.63614691 200000 
261 54431793 261 . 54431024 000000 
43164844506 431.64917860 00016 





be. System 2 


Z é 
as system 
‘h 
—— >= ae mirror image 


System parameters: See system 3, E.21 





Mathematical models: Distributed mass with the 


effects of shear deflection and rotational inertia neglected. 


Results: 





frequency (rad./sec.) 


Difference (7%) 
33044525745 33 44525745 -00000 
108. 38602303 108 . 38602303 »00000 
226.13877715 22613877715 .00000 
386.71069206 386.71059481 -00002 
590.10235871 590.08973387 200211 





c. System 3 7 






System (component 
numbers shown) 
ah 


indicates a pipe 
ca| hanger 


mirror images 


System parameters: 






Component 


H 


diameter (in.) 













wall thickness (in. ) 


length (in.) 


Ba Ba eS 





radius of curvifin. ) 


w 






incl. angle of are 








Intensive properties of all conponents except 
hangers: 
density: 556 lbs./cu.-ft. 
Shear modulus: 12. 10° psi 


elastic modulus: 30: 106 psi 


Hanger spring rates: 
CLK = 1000 1bs7 a 
CLY = CLZ 


LOOO Wes in. 


CTX = CTY CiZ = HCGC-in.-1b./rad. 


il 


End conditions: 


Both ends of system fixed. 


Mathematical model: Distributed mass with the 


effects of shear deflection and rotational inertia considered. 


Results: 
in-plane vibrations 
Mode frequency (rad./sec.) 
116.91660080 116.91660084 00000 
B05ee1 7211521 305.21720649 00000 


B45 .99432421 846 .17624950 02150 
— 113509717500 —— 





1 
2 
3 567 h0373058 567 .56397808 00005 
4 
5 





out-of-plane vibrations 









frequency (rad./sec,) 


system | Mirror Image a. 
tt 
114.21818694 | 114.21818694 



















243 692831220 243 ..92831317 
321 .40891123 32140891770 » 00000 
591 40446121 
872329102583 
196% 01776773 
U5 ol 7 3050 


591 OM O7 74 e 00000 
872 624253187 00131 





1136.196150%5 201567 









in=plane vibrations 


Mode frequency (rad./sec. } 
muiber]  Systen | Mirror Image b. 


al 116.91660080 116,91660080 
e 
3 
4 
frequency (rad./sec. } 


305.21 721131 
Mode 
aunber] System Mirror image b. {Difference 







Difference (/) 








200000 








305221721131 
567 «56373058 
845 99432421 


» 00000 





©90000 





567 056373058 
845 9943421 





» 90000 





out-of-plane vibrations 


1 114,21818694 11421818694 200000 
2 243.92831220 243 .92831220 200000 
3 321 040891123 32140891123 200000 
4 591 40446121 591 40446121 200000 
5 672623102583 O72 se 3102504 200000 
6 1138.01776773 1138.01776773 - 00000 
9 


1571 .78173050 1571 (OL 3050 ©90000 




















d. System 4 


system (component 
numbers shown) 


_Grindicates pipe 
hanger 


mirror image 


System parameters; 









Component 





Parameter 








diameter (in.) 
wall thickness (in.) 


length (in.) 











NZ 3} 
H 
diameter (in.) 250 1c a0 
A 
wall thickness (in.) " Fes ao 
length (in.) A 39.27| 100. 
E 50% | <> 
R 45° 0? 


E-9 








Intensive »roverties of all connonents except 
hangers: 
density: 460 lbs./cu.-ft. 
Shear modulus: 12: 10° psi 


elastic modulus 30 - 10 


Hanger spring rates: 


CLX = CLY = CL 


000. 1bs./ ine 


CTX 


Cig CT2Z 


tOC0> in.e-lb./rads 


Bnd conditions: all nain member and branch ends 


fixed. 


Mathematical model: Distributed mass with the 


effects of Shear deflection and rotational inertia considered. 


Results: 
in-plane vibrations 
Mode frequency (rad/sec.) 


number Mirror Image Difference (j) 


ie ol GeiS2 18.55189208 200016 


pile25211856 51 252127390 -00001 
9094757587 
10719459215 








-_ 
= 
a 





out-of-plane vibrations 


node frequency (rad./sec. 


number Sy ster Mirror image Difference(7) 
14.31754957 14,.31754879 »00000 
34 4 54954K9 34 4 5UONE 54 «00002 
41 ..07914530 41207913499 200002 
7306452245 7306429798 200030 
115.86572509 115.86505834 3 V0057 








E.32 Single Branch System 







ae Member 


Case 1 (component 
nN numbers shown) 


Cp sie 
Main Member 
Case 2 (component 


numbers shown) 
P=Je0. 





diameter (in.) 


wall thickness (in.) 


length (in.) 


madiue Of Curve | in. } 





incl. angle of arc 








Intensive properties of components: 


density: 556 lbs./eu.-ft. 


6 


shear modulus: 12° 10° psi 


elastic modulus: 30> 10° psi 


End conditions: Both ends of main member, and 


end of branch, fixed. 


Mathematical model: Lumped mass with the effects 


of shear deflection and rotational inertia neglected. 


Results: 


out-of=plane vibrations 


Mode frequency (rad./sec.) 
number System 1 System 2 


4964604938 4964604961 
135652684605 135052685066 
269 683567590 269 83569992 
530.72012316 530675718901 
O75 51962199 675057044008 
901 652786818 901 .76873827 























Difference (%) 





»o90000 





000000 





© 00000 









2 00698 












000796 






002671 








In-plane vibrations: 


frequency (rad./sec.) 


132.89829688 














Mode 
number System 2 Difference () 















132 ,89836401 »00005 






190.61015386 190.61009901 200002 
430.03948554 
512 254463160 
893210993095 


96300709091 








430.03806330 
51255370411 
893 «66594383 
962 .07559083 


00033 
.00177 








Ooze! 





009682 





E-13 





4PPeNDIX F 


LRaAwSroR MATRIA SOURCE foro bnCce 


ieee General 

Fourteen of the transfer matrices used in Vier 
were drawn from the catalog of transfer matrices, pages 
376-410, in Pestel's "Matrix Methods in Elastomechanics" 
(a. Some of the transfer matrices constructed in VIPIPE 
are a composite of two sevarate transfer matrices froa[}]. 
These composite matrices are constructed as shown in 
Figs. F.i-1, and F.1-2. The former shows the format of 


case lcomposite transfer matrices, the latter those of Case 2. 





In both 

| j - 

oe : pine Bi 6 

or gal aa: : pee = Agee 
| null 
submatrices 

Oo | A 10 = | 
a 0 A 

—— 7. 


Fig. F.1-1 In-plane composite Fig. F.1-2 Out-of-plane 
transfer watrix format composite transfer natrix 
format 


cases the a matrix is of fourth order, the B matrix of second 
4 


order. In Case ~ composite transfer wsatrices the elenents 


of the B matrix occupy the positions shown (Fig. F.1-1). 





In the following source listing, when two matrices are 
indicated for one subroutine, the first is the A matrix. 
(All component transfer matrices are constructed in sub- 
routines.) 


F.2 In-Plene Vibration Transfer Matrices 





Transfer matrix Subroutine Source 
straight section with distri- DIST C=3a 
buted mass, shear deflection Cola 
and rotational inertia 

straight section, elastic SFIELD C=3¢ 
massless field with shear C-l1c 
deflection 

curved section, elastic mass- CFIszLD C-7c 
less field without shear 

deflection 

point mass with rotational PORE C=-3f 
inertia C=-1b 
rigid straight section with RIGID C=3d 
rotational inertia, distributed C-1b 
nass 

stiff corner (elbow), massless STIFCO C-7f 
field 


F.3 OQut-of=-Plane Vibration Transfer Matrices 


Transfer matrix subroutine Source 
straight section with distri- DISTMO C=-3a 

buted mass, shear deflection C=2a 

and rotational inertia 

straight section, elastic mass- SFIEILO C-3c 

less field with shear deflection C-2c 

curved section, elastic mass- CFILELO C-8c 


less field with shear deflection 


point mass with rotational POINO C-3f 
inertia C-2b 
rigid straight section with RIGIO C-3d 
rotational inertia, distributed C=2b 
mass 

stiff corner (elbow), massless STIFOO C-8f 


frela F-2 


Me Th 
fy | 





APPENDIX G 


SAMPLE PROBL 


G.1 General Remarks 

Appendix G. is intended to supplement apnendix B, 
Included in appendix G. is an illustration of the data 
deck for a typical multiple branch system with hangers, 
and an example of the three ontput formats available in 
VIPIPE. 

G.2 The Data Deck (Tug, G.2-2) 

The data deck as shown calls for determination of both 
Case 1, and Case 2 natural frequencies (10P=0). If only 
Case 1 frequencies were sought, card 4 would be removed 
from the deck, and top! (circled field on the first data 
card) would be set equal to 1. Similarly, if only Case 2 
frequencies were sought, cards 2 and 3 would be removed, 
and IOP set equal to 2. 

In the sample problem configuration all boundary con- 
ditions may be specified by code (B.61), but in order to 
illustrate the method of programming boundary conditions 
for which no code is provided, the Case 1 boundary condition 
state vector for branch 1 is read-in directly. Note that 
the second field of card 2 is blank (blanks are read as 


zeros by the computer) and that card 3 contains the particular 


See appen. B., section B.71. 
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boundary condition state vector. 

The data deck used in Fig. G.2-2 pertains to system 4 
of section E.31, appendix E. The system is illustrated in 
Fic. G.2-1. 

G.3 Output Options 

axamples of the three output formats available in 
VIPIPE, through specification oe the value of the parameter 
GuaBC (B.71), are included on pages G-4, and G-5. The 
system analyzed was system 3 of section 4.31, appendix #. 
and is illustrated there. The problem numbers refer to 
the order in which the data decks are "stacked" (B.8). 

The heading ... "PROGR:aM VIPIPE" ... etc. is printed once 
per run, ¢.g., if 4 data decks were “stacked" behind the 
program deck, and each asked for Case 1 and Case 2 analyses 
the heading would precede the print out of the results of 


the Case 1 analysis of system 1 only. 


G=3 








teheem 


PROGRAW VIPIPE Ge Es FINK WNFA-1 FEBRUARY 1964. 


POCA g N A ANAR PIP 
YETHOO OF TRANSFER ene. EGLENCIES JF PLANAR PIPING SYSTEM ARE ODETERMINEO BY AN [TERATIVE PROCECURE USING THE 


SSSSFCOSSHSSHHHSTH TSA SHCSHSHTCAHHASGHHFHC ESCH CHHTESSHTOHR ETS SBGE SE SOP A CRS HETOTOUD TOAD AR ES SEHD ODSCHAGEAHALASEATDIVOTSEOE 


PROBLEM 1 a 
IN PLANE WOCE FREQUENCIES (RAO ANS PER SECCNO) 
MODE FREGLENCY : 
: 1 116,91660088 
2 305.21720892 
3 0 $67.56629170 
4 . 646.740"5851 


PAOBLEM TERPINATEO CUE TO SICNIFICANT FIGURE LIMITATION OF CCWPUTER. 


A OISTRIBLTEC MASS APPROACH WAS EMPLOYEO ANO 
THE EFFECTS OF SHEAR ofecees ion ANO Rotaay TNERTIA WERE CONS [DEREO. 
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Fig. 63-1 GDABC =O, 
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t  sApGRan VaPEPE : G. Es FIR man] FEBRUARY 1944 
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WETHOO OF TRAMGFER MATRICES. COKE NCHES 2 E PING SY STE Ls : : y ! 


SROSCSSSSSAHSEHSESEDSEFSEFCHASSHSE FHS SGD OSTSE HE SESE CES Sevvesecesepe SOODAHECEES rT SSCSHESHESEHTEHSE SASHES SHSHASHECSHEHSSETCECEES 


s PROBLEM } 

: IM PLANE MODE FREQUENCIES (RADIANS PER SECOND) . 
; MODE FREQUENCY 
' 1 196.91660060 
i F ; 2 | 365521721073 

’ | 5 S67 .56385676 

et h 846 005107199 

| 5 112857638979 
' PROBLEM TERMINATED OUE TO SIGNIFICANT FIGURE LIMITATION OF CCHPUTER. 
| THE EFFECTS OF SHEAR DEFLESTION ANG ROTERY. INER TA WERE CONS IOERED. - 


Ee" tial aaa aa Oa lahat EI I aaa i 
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-000 +0 139.900 3 200 a) o3CE408 - . 
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°00 b 398 3 200 Joes 12669 
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' § £000 : £388 -b3000 2 * 90. 25:88 7356288 15eeda 
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3 & _ 1000.00 1000.00 1000 00 1000.00 1000.00 1000.c0 
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‘ wie es e @ |: {: . 
SVE= » =| © Vo Be Ie 


mueveNee? CDABC= 1. 


6-4 





2 


PRORLCM 
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“FREQUENCY 


MODE 
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INPUT DATA 
COMPONENT 


LENGTH 


WALL TRICKNESS 


OLAMETER 


Stet i "La 


HANGERS 


cCLY cur cix cry cr2 
1000.00 1000.00 


cLx 
1000.0¢ 


COMPONENT 


yoc0.co 


1000.00 
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